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PRESENTACION 





Esto ejemplar de la Revista SELPER forma parte de 
una serio de cuatro (4) números que han sido lanza- 
dos, con jerzo, en un Plan Editorial de 
Emergencia, destinado a recuperar la continuidad 
editorial de nuestra Sociedad, iniciada durante la 
¡Sede deSELPER en Chile (1983-86), la que, lamenta- 
blemente por motivos de fuerza mayor, no se pudo 
seguir en la forma en que hubiésemos deseado, Es 
así como se han lanzado números especiales fecha- 
dos simbólicamente como Sealembre 1960, Sep: 
tiembre 1987, Septiembre 1988 y Septiembre 1989, 
para representar la continuidad en esos años y co- 
menzar a publicar normalmente en 1990. Dado que 
¡son Revistas editadas con posterioridad (Septiem! 
1989, facha do matorialización cel Plan Editorial de 
han incluido Secciones de Noti. 
Iberoamericanas e Internaciona- 
aparecerán en los números regulares 
blicados durante 1990 en forma semes- 
iberoamericano e internacional es 
adecuado, en forma trimestral 

Los artículos técnicos seleccionados para la Rovi 
ta SELPER 1986, corresponden a trabajos preser 
dos en ol Seminario de Geología Regional co-organi- 
zado porla Unión Internacional 
cas (1UGS), UNESCO, SELPER 
quisas Espacial (NDE) de Bra (Sao José, dos 

'ampos, SP, Brasil, Diciembre 1985). La continua: 
ción de estos anículos técnicos conforman la pane 
medular de la Revista SELPER 1987, Par 
SELPER 1988 se ha seleccionado un ¡me y 
completo artículo del Dr. Jorge Lira (Universidad 
¡Nacional Autónoma de México, UNAM), sobre Carao- 
torización Matemática de Imágenes Digitales, el que 
sehacraido conveniente complementar con su: 
mido pero ejemplficador trabajo del Ing. Marcelo 
'Campl (Comisión Nacional de Investigaciones Espa- 
ciales de Argenina, CNIE), destinado a profesionales 
o tan especializados, Finalmente, para la Revista 
SELPER 1989, se continuó con el siguiente capítulo 
escrito por el Dr, Jorge Lira (UNAM), sobre Transfor- 
maciones Reversibles y Realces Globales, Espera- 
mos que estostemas sean de nterós y recibira futuro, 
para las próximas Revistas SELPER, otras valiosas 
contribuciones de la comunidad profesional latinoar 
Americana y mundial, 

Probablemente llame la atención el por qué de este 
sfuorzo yla oportunidad en que éste so materializa, 
Las razones son profundas y constituyen un compro- 
miso dehonor, decontinuidad y constante perfeccio- 

ctividades de SELPER alolargo de 
lo de SELPER 
sn Brasil (1986-89), se obtuvieron inmensos logros.en 
amplar la cooperación intemaciona, organización 
de Simposios Latinoamericanos en diversos países 
de Nuestra América en coordinación con la Sede 
SELPER en Brasil, reconocimientos intemacionales 
hacia nuestra Sociedad y diversas actividades de 
apoyo dentro y fuera de Iberoamérica. Todo ello 
demuestra la madurez y consolidación que ha logra- 
doSELPER. Lamentablemente, por razones de fuer. 
za mayor, no fue posible reforzar la actividad eciorial 
de SELPER. 

Conscientes de la importancia que tienen las Publi 
caciones como agentes de comunicación y perfec- 
cionamiento permanente de nuestros Socios, ya que 
no es posiblo para todos participar en los diversos 
Cursos de Perfeccionamiento, Seminarios, Simpo- 




















































































sios y otras actividades organizadas o co-patrocina: 
porSELPER y, conia misma visión y generosidad 
con que la Sede SELPER Brasil abordó la realización 
de Reuniones Plenarias y Simposios Latinoamerica: 
nos fuera de su país, también ha.existido una apertura 
para realizar un Plan Editorial a más largo plazo, 
solictando y traspasando importantes responsabil:- 
dades al Comité Editorial que tantos éxitos logró 
durante la Sede de SELPER en Chile. Ello, sumado a 
tin representativo Comité de Expertos y de Arbitraje a 
nivellatinoamericano, compuestos por profesionales 
de probada excelencia y dedicación a SELPER, junto 
ala esperada decisión de las futuras Sedes de conti 
nuar con un Comitó Editorial más estable -indepen- 
dientemento do la rotación de óstas. deberá pormitir 
lograr la tan ansiada continuidad editorial que tanto 
necesita nuestra Sociedad para fortalecer y multiplk 
car aún más los éxitos cosechados en todos estos 
años, desde sus inicios en Ecuador, pasando por 
Chilo, Brasil, Argentina y posteriormente otros palses 

hermanos. 
Deseamos agradecer la valiosa colaboración presta- 
da por los diversos Socios e Instituciones de los 
Capitulos Nacionales SELPER, en particular los de 
Argentina, Bras, Chile, Colombia, Ecuador, México, 
Pera y ojos ptes hermanos de América del Sur y 
rica Central, sin cuyo apoyo y participación per- 
mmanonto, SELPER no hubleso sido. una realidad. 
los Socios, Instituciones y Empresas de 

Honorarios de Alemania, Canadá, 
YU, Francia, Holanda y Portugal, Igual: 
), las organizaciones internacionales como la 
Agencia Escacial Europea (ESA), Naciones Unidas 
(ÓNU), IUGS, COSPAR/ICSU y otras, que tan decisi- 
vamente también han aportado a los logros de SEL- 
PER, Finalmente, para la realización de este Plan 
Editorial de Emergencia, vayan nuestros sinceros 
reconocimientos al Comité de Publicaciones SEL 
PER Chile, Editorial Universitaria, Impresora América 
y, Por supuesto, al Insttuto Geográlico Miltar (GM) 
de Chile -en especial a su Director, Brigadier Don 
Roberto Castillo Vie, subdirector Técnico; TCL Lauta- 
roRivas, TCL Sergio Matus y MY. Mario Molina, pues 
sa colaboración será vital para consolidar este Plan 
Editorial a largo plazo de SELPER. 

Do esta forma, os un gran orgullo para la Sede 
'SELPER Brasil poder entregar una situación Editorial 
aldía altérmino de su gestión; para SELPER Chi 
¡gran honor poder continuar contribuyendo en forma 
lan decisiva a la mayor consolidación de SELPER; y, 
para SELPER Argentina, una gran satisfacción y un 
desafío el poder continuar avanzando caca vez más 
mn estas actividades en beneficio de SELPER y sus 
ideales de cooperación y hermandad iberoamericana 
y universal. 
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EDITORIAL 


Thisissue,thatis parto! a series offour(4) SELPER 
Roviaws, constítutes a great effor to continue withthe 
Editorial Programme that was begun by SELPER- 
Chile during 1983-86. Unfortunately and due to seri- 
0us problems, this could not be continued. Because 
otthis, we are publishing now these special editions 
which will represent the continuy along years 1986- 
1989. These issues have been symbolically dated 
September 1986, September 1987, September 1988 
and September 1989. Due to the fact that these Re- 
viaws were edited in September 1989 (date of the 
materialization of this Emergency Editorial Plan), ner- 
thor international and Iberoamerican Activities. nor 
News were included, These Activities and News are 
ing to appear in the normal sections in the regular 

which will be published during 1980, every six 
months but if the Iberomerican and Intemational 
Supporia utablo, SELPER Review published 
Quarter 


“To technical papers selected tor SELPER Review 


1986, come from the Seminar on Regional Geology, 


Eoprganized by Ineratonal Union ol Geological 
Sciences (UGS), UNESCO, SELPER and INPE 
(Brasilian Space Research Institute), held in Sao Jose 
dos Campos, Brazil, on December 1985. SELPER 
Review 1987 was shaped with most ofthese papers 
avallablo at tho Editorial Committee. An important 
paper about Mathematical Charactorization of Digital 
Images by Dr, Jorge Lira (UNAM, National University 
ofMoxico) makes SELPER Review 1988 togetherwih 
a summarizod paper by Eng. Marcelo Campi (CNIE, 
Argentinian Space Research Commission). Finally, 
SELPER Review 1989 contains the next chapter w- 
ton by Or. Jorge Lira (UNAM) about Roversible Trans- 
formations and Global Enhancementes, We hope 
those papers will be of interest and hope to receive 
morecontributions forthe next SELPER Reviews from 
Latin American and worldwide community in the next 
future 

The reasons to do this Emergency Editorial Plan, at 
this timo, constitute a commitment and a constant 
improvement ín the activities of SELPER all along the 
rotation of ts Seats. In Brazil (during 1986-69), 
SELPER got alot of achievementsinthe international 
opening and coorganization ot Latin-American Sym. 
posia «which took place ín different Latin American 
countries. Our Society also received additional iter- 
national recognition, as wellas international support, 
The above showed that our Society had reached an 
important consolidation and maturity. ln spite ofthis 
was, Unfortunately, not possible for our Society to 
reinforce SELPER's Editing activities. 

Aware of the importance that these publications 
haveas agents ofcommunication and constanttoolof 
improvementfor our members, sinceitisnot possible 
for all of them to take part in the different training 
Courses, seminars, symposia or anyother activites 
which are organized or co-sponsored by SELPER 
and, with the same view and generosty which 
SELPER/Brazll_ organized Plenary Meetings and 
Latin-American Symposia with abroad, there has also 
existed 1he same wil to accept and to create a long- 
term Edition Program by applying and transfering 
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important and significant responsibilties to the Edit- 
ing Committee that had so successful participation at 
he SELPER Seat in Chil (1989-86). Althis, plus the 

'ofboth a Specialists Committes and Rete- 
ree! thatare consítutedby professionals ol 
high excellence and devotion to SELPER, together 
wi the expected decision of the future SELPER 
Seatsinorderto continue workingwih amore perma- 
ent Editing Committes thatis not going to be affected 
bySELPER Seats rotation, we hope that twil allowto 
getthe so expected editing continuiy that our socioty 
eds to reinforcs and muliply even more the 
achievements got in the past, since is begining in 
Ecuador, going 1hrough Chile, Brazil, Argentina and 
other brother countries in the nar future. 

We also want to thank the valuable cooperation 
given by the difterent members of SELPER and the 
Iremos ver 1 National Chapter ol SELPER 
worked specially those from Argentina, Brazil, Chile, 
Colombia, Ecuador, México, Perú and other brother 
countries of Souh' America and Central America, 
¡whose support and cooperation and permanent par. 
ticipation were fundamental forte creation and exis- 
tanos ol SELPER, We aso vt 10 thank to alte 
Members, Enterprises and Institutions of 
Honorary' Chapters, special those of Canada, 
France, Germany, Holland, Portugal, Spaln and 
USA Ao o ose earaienal oraciziora 

'uropean Space Agency (ESA), United Nations, 
1UGS, COSPAR/ICSÚ and others who have been 
decisivo for SELPER success. Last but no least, wo 
want to thank the valuablo support given by tha 
'SELPER-Chile Editing Comminee, Editoral Universi- 
taria, Impresora América and Instituto Geográfico 
Miitar (1GM) from Chile - special to ts Director, 
Brigadier General Mr. Roberto Castillo V.; Technical 
Subdirector, LTC Lautaro Rivas; LTC Sergio Matus 
and MY, Mario Molina whose colaboration will bo of 
real importance to pertorm this SELPER's long+orm 
Edition Programme. 

pi ey Ll al prou for SELPER Era] to 
offer an updated edition programm at 1he end of ts 
mandato;  ís a real proud for SELPER-Chile to con- 
tínue: nsuch a decisivewaytothe greater 
Consolidation ot SELPER; is a real challenge and 
satisfaction for SELPER-Argentina to continue im- 
provingthese activities for the benefit of SELPER and 
s important aims of Iberoamerican and Universal 
brotherhood and cooperation. 






(1986-1989) 


Prot. MAURICIO ARAYA-FIGUEROA 
'Vice-President (Ecuador, 1980-83) 
TES arca 
SELPER Seat Brazil (1986-89) 

)- MIGUEL SANCHEZ- 
posa erro 

EN 

Sl 
ona a 





SELPER septiembre 1989 





TRANFORMACIONES 
REVERSIBLES: REALCES GLOBALES 





Dr. Jorge Lira 


Instituto de Geofísica. Universidad Nacional 
Autónoma de México (UNAM). A. P. 22-560. 
Deleg. Tlalpan. 1.400 México D.F. MEXICO. 


Resumen 


Eneste completo trabajo de revisión se continúa con eltema iniciado en 
laRevista SELPER anterior, constituyendo esta vezel capítulo 6 dellibro 
“Introducción al Procesamiento Digital de Imágenes” del mismo autor, 
razón por la cual se ha mantenido la numeración original. Para una 
mejor exposición, el tema se ha subdividido en: 6.1) Consideraciones 
Generales; 6.2) Transformaciones en el Espacio de Coordenadas; 6.3) 
Modificaciones al Histograma; 6.4) Operadores Diferenciales; 6.5) 
Transformada de Fourier; 6.6) Descomposición en Componentes 
Principales; 6.7) Análisis Canónico; 6.8) Transformación Kauth- 
Thomas; 6.9) Transformaciones Geométricas; 6.10) Transformaciones 


Quasi-lineales; Referencias. 











| realce global de la imagen multiespec- 
E tral, banda por banda, o en su conjunto, 

tiene como objetivo principal el de facil 
tar la apreciación general de la misma 








"por medios visuales y digitales. Una imagen global- 
'menterealzada puede ser empleada comotal para un 
reconocimiento analógico de patrones, o bien vista 
¡como un paso previo para un reconocimiento digital 


6.1. Consideraciones Generales 


Enslcapítulo anteriordejamos entendida laidea de 
imagen multespectral, los elementos que la compo- 
nen, lainterrelación que hay entre ellos yla estructura 
¿que forman. Muy vasto es el campo que se refere aa 
estructura de una imagen y para atacario tendriamos 
que entrar de lleno en el terreno de la psicovisión, 
“amén de que mucho de esto se refers al reconoci- 
miento de patrones, lo que está fuera del alcance de 
esta obra, Aun así sabemos ya que es una imagen 


digital mutiespectral y podemos proceder ahora a 
estudiar qué podemos hacer con ella, es decir cómo 
la podemos transformar para extraer la información 
perinente a nuestros objetivos. El primer gran grupo 
de transformaciones reversibles, incluyendo las uni- 
tarias y que por su misma naturaleza permiten recu- 
prar,en principio, laimagen original. En lastransfor- 
'maciones reversibles existe siempre el operador 
verso que permite recobrar la imagen original a partir 
de la imagen transformada. Como ya se mencionó 
anteriormente, a estas transformaciones se les cono- 
ce también como realces generales o realces no- 
selectivos puesto que no se realza una clase de 
patronesacosta deotra. Esto auncuando enalgunos 
casos ciertos patrones sean más realzados que 
tros. Algunos operadores como los de tipo diferen- 
cialtiendenarealzar con mayor énfasisciertosrasgos 
de la imagen, pero esto no es a costa de los otros 
rasgos y menos aun ningún aspecto de laimagen es 
suprimido. Por tanto, incluyendo estos casos, esta- 
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mos hablando de realces no selectivos ya que no se 
está segmentando selectivamente la imagen. En un 
(salce general, al no segmentar la imagen, sino pre- 
servar todos los patrones que la componen, reaiza- 
dos en menor o mayor grado, se genera una imagen 
demejor calidad para propósitos específicos. Deesta 
manerase conserva todala información de laimagen 
original y por esto setiene una transformación rever- 
sible, Podemos deci, en otras palabras, que desea- 
mos realces globales para hacer una apreciación 
¡general de la imagen, en una primera etapa de anái 
sis y esto implica que los realces sean no-selectivos. 
La principal ventaja de que una transformación sea 
reversible es que podemos establecer otra represen 
tación de laimagen con el objeto de visualizar ciertas 
caractoríticas o estructuras de la misma. Una vez 
visualizadas estas estructuras, podemos regresar ala 
imagen original (si así conviene) y establecer útles 
conclusiones, Una segunda ventaja de este tipo de 
roalces generales es que son de rápida ejecución 
desde el punto de vista computacional. Además tie- 
¡non pocos o ningún parámetro de entrada; todo esto 
las hace de fácil uso y sencilla implantación en com» 
putadora. Algunos do losrealces generales se aplican 
bandaporbandayotros a laimagen multiespectralen 
su conjunto, Las modificaciones al histograma o la 
transformada de Fourier son ejemplos de los prime- 
109 yla descomposición en componentes principales 
ola transformada de Kauth y Thomas son ejemplos 
delos segundos. Al realizar transformaciones banda 
por banda, tenemos la opción de varía el valor de los 
Parámetros (si es que existen) involucrados realzan 
dolas bandas en diferente grado. Como cada banda 
roprosenta diferentes aspectos de los patrones quela 
componen, al aplicar el realce en cuestión se entati- 
'an entonces diferentes rasgos de la imagen por 
separado, Entransformaciones aplicadas ala imagen 
multiespectral como tal, el realce se realiza a part de 
todas las características presentes en cada banda, de 
staformaciertas clases de cobertura pueden quedar 
evidenciadas como tales, Esto os ol caso de la des- 
composición en componentes principales y de la 
transformada de Kauth y Thomas. Para ciertos pro- 
blemas los realces banda por banda son adecuados, 
paraotros problemas los realces de conjunto son los 
mejores, Esto será discutido con mayor amplitud en 
las secciones correspondientes a cada transtorma- 
ción. 

Existen varias formas de ver a una transformación 
reversible. Laprimeraes quelatransformación de una 
imagen puede interpretarse como una descomposi- 
ción de ésta en un espectro generalizado bidimen- 
sional empleando una base de funciones ortonorma- 
los, Cada componente espectral en el dominio de la 
transformada corresponde a la cantidad de energía 
de la función espectral dentro de la imagen original. 
En este contexto, el concepto de frecuencia puede 
ser generalizado para incluir transformaciones de 
funciones diferentes alas de Cosenoy Seno. Estetipo 
de análisis espectral generalizado es útil n la inves- 
tigación de sintetización de estructuras específicas 
para una clase dada de imágenes. Orra forma de 
visualizarla transformación de una imagen, es consi- 




















derar ala transformación como una rotación mutidi- 
mmensional de coordenadas, es deci, como un cam- 
blo de la base que genera el espacio al cual se 
“encuentra referida la imagen muttiespectral. Una ter- 
¡cera manera de ver esto, es considerar a la transfor- 
mación como un medio de sintetizar una imagen con 
una función matemática brdimensional para un 
dominio fo de transformación de coordenadas, En 
esta interpretación, el operador de la transformación 
ss llamado función base bidimensional y elcoeficien- 
to delatransformación, (4,1), » (4,1), esla amplitud 
delas funciones baserequeridas para la síntesis dela 
imagen. Cabe mencionar aquí que una de las princi- 
pales propiedades de una transformación reversible 
s que se preservan los aspectos topológicos de la 
imagen original, En algunas transformaciones con 
mayor grado de simetría (unitarias) se conserva, 
además de lo anterior, la métrica. Ahora podemos 
¡decir que, en términos generales, el principal objetivo 
de unrealce es el de presentar una imagen multas. 
pectralenunmarco de referencia más adecuado para 
una aplicación específica, lo que quiere decir que el 
tipo de realce está definido por el problema en cu 
tión, Es así que unaciertatécnica de realce puedo ser 
adecuada para una radiografía, pero no para una 
¡imagen sateltaria. No existe unateoría ger 
realce de imágenes, pero sí puede afirmarse que el 
diseño de un realcetiene como propós 
tal el de poner en evidencia uno o vas 
espaciales o espectrales, Lo que es subjetivo es la 
evaluación de qué tan bien un cierto realce haevide 
cado un determinado patrón, Para minimizar esta 
subjetividad, es necesario predefinir un marco de 
fidelidad al cual se refiera la imagen realzada, aun así 
no existen criterios generales para establecer marcos 
de ficelidad, sino que más bien éstos están en función 
de una aplicación específica y del valor que se le 
asigne a la imagen realizada. Ahora podemos decir 
Quelos reales se dividen en dos grandes categorías: 
aquellos que se basan en la expansión de funciones 
ontonormales como la transformada de Fourier, de 
Hadamard o de Haar. La otra categoría es la que se 
roñere ala ejecución, en el espacio de coordenadas, 
dela manipulación directa de los pixels de la imagen, 
ya sea a través de ventanas o de los pixels de la 
imagen en su conjunto, Como ejemplo de realces 
llevados a cabo en el espacio de coordenadas tene- 
'moslas elongaciones al histograma y los operadores 
diferenciales aproximados por diferencias finitas. En 
resumen puede afirmarse que los realcos, análog; 
mente a tras transformaciones a la imagen, puede 
ejecutarse en el espacio de coordenadas o 
frecuencias generalizadas. 
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Transformaciones en el Espacio de Coor- 





Deacuerdo alo visto en elcapitulo, el proceso de 
digitización de una escena produce un arreglo de 
números en forma matricial. Seal (1.1) la matriz que 
representa una banda cualquiera de una imagen 
digital mutiespectral. 
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es decir, que esta imagen tiene N líneas de M pixels 
cada una. Recordemos que cada elemento de esta 
matriz representa univocamente las características 
promedio dela radiación proveniente de los CIVde la. 
escena, para un cierto intervalo de energía. Una 
transformación general discreta de la matriz. (1-11, 
esuna matriz. (+) de la misma dimensión y definida 
¡como el producto de tres matrices. 


(FI teIteItO) 


dondelas matrices. PI y 10). sonnosingulares* 
y d dimensiones (M,M) y (N.N) respectivamente. La 
transformación de la imagen ¡16.311 queda 
determinada una vez que las matrices (+! y (01son 
definidas. Cómo establecer definiciones apropiadas 
para. (?1 y 19] esuntema extenso que puede ser 
estudiado (Arfken, 1970) en la Iteratura. Cabe ahora 
mencionar que el hecho de que (7) ; 10) seanno- 
singulares hace que la transformación ariba definida 
(ec.6.2) seainverible y puede puntualizarse también 
que sino existen los inversos respectivos a (PI y 
101 se pueden estimar sus pseudoinversos (Pratt, 
1978). En otras palabras, aquellas transformaciones 
lineales y no-nvertibles, implican entre otras cosas 
quelas matrices (Py (0) - sonsingulares, aunque 
si bien esto es una condición necesaria más no 
suficiente, puesto que — [FI y 10) pueden ser 
singulares, pero bajo ciertas restricciones es posible 
encontrar su pseudoinverso. En forma explícita, la 
transformación discreta de *(:.)1se puede escribir 
como sigue 
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parak,1 = 1,2,..MN respectivamente. Silas matrices 
(PI y 101 de latransformación sonno-singula- 
res, entonces sus Ínversos están univocamente 


definidos. Sean (?)-' + 101'losinversos del?] y 191 
respectivamente, entonces multiplicando a: en (8.2) 
porlaizquierda por '7¡-' y porla derecha por (0), 
se tiene 


que es justamente la transformación inversa, es decir 
:18s la transformación inversa de ''| . Una manera 
"conveniente de escribir la forma explícita de la trans- 
formación discreta es la siguiente 








y apartir de aquí se define lo que se conoce como 
kernel de la transformación 


es decir, que!1*! es una matriz de dimensión (M.N) 
cuyo elemento (m,n) es precisamente I?k.s0(0.1). 
Elkemel B de la transformación inversa puede obte- 
nerse en forma similar. En una transformación unita- 
fía, los kereis A y B cumplen con un juego de 
condiciones de ortonormalidad. Este juego es el 
siguiente 














donde elsimbolo * significa conjugado complejo yla 

* serefere ala deta de Kronecker, En estas con- 
diciones está claro que los kernels A y B deben incluir 
el factor de normalización correspondiente. De 
regreso alas relaciones (5.5) y (6.6) podemos escribir 
la transformación a una imagen como 





« Presentamos equí on forma resumida las propíecados matriciales empleadas an esta sección: a)- Una matiz cuacrada posee 
pseclolnwers ss determinante es diererte de Cero, als marcas se aman no-sinutres; 1) Una matt real, cuadrada T os simélrica 
EITOZ Y E donde T.* esla tanspuesta de Y; c)- 1 esotogorsísl 7.1 7 = /, donde | esta matriz idontdad, de aquíque 
ara una mabtz ro uncrad que 308 smméerca y ortogonal que ene. T =1 = T, d)- Una mare compa cuactada O, es Hlermina sl 


(3 = € donde O "es el conlugado complejo de € 
que oa Htrmilana y unitaria se eno CF = O 





€ esuntada sl CC = 1 ¿do oquíquepara una matiz cuadrada compleja 
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ques finalmente un modelo matemático detranstor- 
mación de imágenes, similar al desarrollado en el 
capítulos (ec.5.50) pero concebido en forma dieren- 
te, En el contexto actual la relación (6.8) puede verse 
de la siguiente manera. 





es decir, que el kemel A es un operador actuando 
sobre la imagen de entrada para producir la imagen 
¿de salida F. Aun cuando la relación (6.8) es básica: 
10 la misma que la (5.50) el contexto bajo el cual 
o aplican es diferente, La relación (5.50) es general y 
¡entiendo en un contexto amplio, no así la relación 
(6.8) que normalmente se aplica cuando la transtor- 
mación se hace en el espacio de la imagen. La 
Coctumbre ha hecho que al término kernel se refiera 
aoperadores que actúan en el espacio de coordena- 
das. Eloperador kernel puedo ser de respuesta finita, 
sobre una ventana, como en el caso de un operador 
diferencial aproximado por diferencias finitas, o pue- 
¿de sor de respuesta infiita (Cappelini, 1978) sobre 
todo el dominio de la imagen. En algunos casos se 
confunde (erróneamente) el kernel con la ventana 
misma. La ventana oslarogión sobre lacual se aplica 
el kernel, es decir éste es la regía bajo la cual se 
combinan los valores de los pixels de la ventana para 
substituir a uno de ellos determinado de antemano. 
En forma simbólica esto se escribe como 











donde (4,1) son las coordenadas del pixel predeter- 
minado, son los pixels de entrada, ps es el pixel o. 
valor de sálida, K es el operador kemel y V., es la 
ventana de dimensión (m n), que puede inclu? desde 
unsólo pixelhasta la imagen entera. Sik significa una 
operación de promedio sobre la ventana V,., enton- 
ces latransformación no es reversible, Por Stro lado 
debemos decir que la ventana tiene que ser una 
región conectada simple, en principio de forma arbi. 
traria y que contiene un cierto número de pixels de la 
Imagen original. Es posible también concebir hiper- 
ventanas definidas sobre la imagen mutiespectal. 
Estas consideraciones implican que los índices m y n 
de la relación (6.8) no necesariamente van hasta la 
dimensión (MN) de la imagen. Veamos ahora que, 
cuando hablamos de un kernel pensamos desde el 
punto operacional, mientras que cuando nos refer- 
mos a la función de respuesta al impulso queremos 
dar una connotación más bien física aa transtorma- 
cióndeunaimagen. La consideración de unaventana 


finita, es decir menor a la dimensión de la imagen, 
implicauna aproximación alatransformación respoo- 
tiva, esto le quita en términos estrictos pero reales, 
reversiblidad a la misma. Hay ocasiones en las quel 
kernel de la transformación no es evidente o no tiene 
forma explícita. En otras ocasiones se parte de una 
forma explícita y debido a ciertas aproximaciones 
introducidas por conveniencia, llega a perderse la 
videncia del operador en la forma dada por la ecua- 
ción (68). Es conveniente ahondar en qué se quiere 
decir con realces en el espacio de coordenadas o en 
el de frecuencias. Un realce en el espacio de la 
imagen implica que la transformación realiza un 
mapeo del espacio de coordenadas sobre sí mismo, 
sn todo caso cambia la escala y la distribución de 
valores de los pixels de la imagen, pero la represen 
tación de la misma no cambia. En los realces que se 
reñeren a espacios de frecuencias generalizadas la 
representación de la imagen sí cambia, En este caso 
la imagen es enviada (homeomorfismo) de un espa- 
¡slo a otro dado que en éste se espera realzar o poner 
¡en evidencia ciertos aspectos que no pueden obser- 
varse directamente en el espacio original, Así, por 
jemplo, la transformada de Fourier de una imagen 
Puede mostrar ciertas propiedades geomélricas dela 
misma, tales como aspectos de simetría de su estruo- 
tura. La distribución de frecuencias de los valores de 
los pixels de una imagen constituya un aspecto de su 
queno puede apreciarse en el dominio de 
sino sólo por medio de la transformada de 
Fourier. Otros aspectos geométricos podrán visual 
zarse con diversas transformaciones do frecuencias 
generalizadas. Como puede verso, un realce general 
serealiza con el objeto de poner en evidencia ciertos 
patrones, o clase de patrones, o la morfología de 
stos, sin detrimento de otros y llevada a cabo esta 
operación, ya sea en el espacio de coordenadas o an 
el de frecuencias generalizadas. 


6.3.- Modificaciones al Histograma 


Uno de los realces más sencillos que pueden apli- 
carse a una imagen es el que se refiere a la modíf- 
cación de su '0 función de densidad. El 
principal efecto de asta transformación es el realce de 
contraste de la imagen, debido a esto es necesario 
revisar primero qué se entiende por constraste en una 
imagen digtal. En cualquier caso el contraste se 
rafiere siempre ala relación de brillantez do un objeto 
¡con su entomo, es decir que de cierta manera se 
relacionael valor de un pixel o una clase de pixels con 
los vecinos próximos. En una apreciación cualitativa 
pormecios visuales, el contraste puede depender de 
las condiciones especificas de iluminación y de facto- 
res particulares en la respuesta del sistema opto- 
neuronal del ser humano. Al definir factores de con- 


an os factores de contraste y siesta modificación es 
'en el sentido apropiado, decimos entonces que el 
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Fig. 6.1. Modilicación de los valores de 
entrada o de acuerdo a la regla T. 





contraste de una imagen digital ha mejorado. Vea- 
mos ahora, entonces, con detalle las implicaciones 
concretas que tiene la modificación del histograma y 
cuales sonen consecuencia los crterios para estable- 
eras transformaciones más adecuadas. Ya vimos 
nel capítulo 5 que ol histograma, bajo ciertas condi. 
ciones estadísticas, es una buena aproximación de la 
unción de densidad (ec. 5.68). Al modificar el histo- 
grama, modificamos la función de densidad y cambia 
en consecuencia la distribución de valores de los 
pixels de la imagen y esto implica necesariamente la 
variación de los factores de contraste de la imagen 
digital. Veamos entonces la metodología general 
para modificar el histograma de una imagen digital, 
paraesto sean. y - los niveles delos pixels de las 
imágenes original y realzada respectivamente y su- 








pongamos que para cada pixel - en la imagen 
original creamos un pixel - enla imagen realzada de 
acuerdo a la siguiente regía 


un ejemplo de cuya forma está dado en la figura 6.1. 
Elefecto de estatransformación para el ejemplo de la 
figura 6.1, es el de producir una imagen de mayor 
«contraste en la parte intermedia de los valores de los 
pixels y que hace más oscuros losniveles bajos y más 
claros los nivelesaltos. Los niveles que se encuentran 
pordebajo del valorm son oscurecidos, mientras que 
aquellos que se encuentran por arriba de m son 
abrilantados. Este tipo de realcs que implica desde 
luego una modificación al histograma, se conoce 
¡como elongación de contraste o elongación del his- 
tograma. Visto de otra manera este realce, podemos 
decir que tanto los niveles bajos (oscuros), como los 
altos (claros), son comprimidos en un rango menor 
que el original. A medida que esta compresión se 

hace más acentuada, detal forma que en ellímite esta 
transformación (fg.6.1) produce unaimagen binaria, 








de manera que se produce una segmentación de la 
imagen cial es deci que esa case de segmen 
tacón es un caso paniular de una elongación de 


De las técnicas de realce más soconídas por su fácil 
valuación computacional y por os resultados visua- 
les que proporciona, es el quese rañere a la modilica- 
ción lineal o no-ineal del histograma de valores de 
una imagen digital; es decir, el histograma es 
selongado lineal o no linealmente. En el caso de una 
imagen mutiespectral, la elongación serealiza banda 
por banda y puede ser diferente para cada una de 
ellas, aunque podría considerarse una 
multielongación que tendría en conjunto los mismos 
resultados que por separado cada banda. El 
histograma representa desde el punto de vista 
estadístico, una propiedad promedio de la imagen y 
como tal proporciona una descripción global de la 
apariencia visual de la imagen. De aquí que la 
modificación adecuada del histograma puedo 
cambiar dicha apariencia con benelicios de realce 
para un propósto específico. Regresemos ahora a 
nuestra” representación matemática de la 
modificación al histograma o elongación de contras- 
to. Para esto supongamos que el rango permitido de 
valores, tanto para la imagen original como para la 
realzada, es el mismo 





sto para todas las bandas de las respectivas imágo- 
nes mutiespectrales. En la práctica, vi» 1 y vf-2n, 
con n = 7 u 8. Se supone, desde luego, que la 
transformación inversa existe y se denotacomo sigue 
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Silatranstormación es lineal, la función T es univa: 
luada y monotónica, sila transformación es no-Íneal 
entonces T es claramente multivaluada. Cabe aclarar 
'quenosevaa considerar cualquier relaciónno-ineal, 
sino únicamente las trigonométricas puesto que és- 
tas tienen bien definido su inverso. Por otro lado o 
y 1 son variables aleatorias discretas dofinidas en el 
intervalo! »..»,)ya que representan los valores de los 
pixels de la imagen. De esta forma, la distribución de 
valores de o y 1 puedeser caracterizada atravós 
de las funciones de densidad 7a(9) y pto 
respectivamente. Es apanirde estas funciones, como. 
yasemencionó, que se pueden obtener propiectades. 
de conjunto de la imagen, en este caso concreto se 
pueden obtener conclusiones acerca de la 
distribuciónde valores de los pixels de la imagen. Por 
ejemplo, enlafigura6.2ase muestra el histograma de 
una imagen que no tiene predominancia de ningún 
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una densidad tipo Gaussiana. 





Y mM y 


Fig. 6:2. En la izquierda so muestra un ejemplo de histograma con densidad uniforme y en la derecha 








tonos oscuros y claros, con una predominancia entos 
tonos intermedios. En otras palabras, a partir de un 
histograma, podemos deducirsila imagentiene poco 
istograma estrecho), buen contraste (fig. 
6:28), o si la imagen tiene o no predominancia de 
cierto conjunto de tonos. Veamos ahora como 
podemos relacionar el de la 
original con aquel de la imagen realzada a través de 
laregla detransformación T. Para esto consideremos 
la Kg. 6.3, en donde se aprecia al histograma de la 
imagen Original y por medio de la regla de 
transformación Y construyamos gráficamente el 
histograma de la imagen realzada. Desde luego, lo 
importante aquí esla selección dela regla para que 
58 produzca el realce que convenga a una aplicación 
específica. Para expresar matemáticamente loque se 
representa gráficamente en la fig.6.3, consideremos 
Que el área bajo la curva de los histogramas original 
y modificado debe ser la misma por las razones de 
conservación mostradas en latabla 5.1, de tal forma 
que se tiene la siguiente relación 





50 
de tal manera que en el límite, cuando ss y 5.0 
Se tiene. 
pao roo 


Ahora bien, recordando que r; T(0) y que o» T! () 
y substituyendo esto en (6.14) se llega finalmente a 
una relación ampliamente conocida en cálculo dife- 
rencia 





alo) 


O tr 


(6.16) 


y que constituye finalmente una expresión analítica 
para el histograma de la imagen realzada; desde 
luego se requiere que exista T-! y esta es una razón 
más por la cual consideramos únicamente transfor- 
maciones reversibles. Es posible que en algunos 
casos no sea factible obtener T-! en forma anallica 
pero sí en forma de tabla; de una u otra manera 
ónemos que suponer que existe este inverso. Otra 
forma conveniente de escribir a relación (6.16) es la 
siguiente 
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esta Última expresión es una relación para transtor- 
mación del histograma en general. Veamos ahora 
varios casos específicos detransformación anotando 
que el usuario puede fácilmente diseñar la propia. El 
Primer caso, el más sencillo, es aquel que se reñere a 
la elongación lineal del histograma o realce líneal de 
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¡contraste como también le Nlaman algunos autores. 
En muchas situaciones, debido a la geometría part- 
cular: fuente/objeto/sensor y a la respuesta particu- 
lar del sensor, no se explota todo el rango dinámico 
de valores disponibles, es decir los valores de los 
pixels no cubren enteramente el intervalo: O -2?- 1 y 
esto hace que la imagen no tenga todo el contraste 
que podría tener, en estos casos es adecuado em- 
plear una elongación lineal del histograma. Para este 
tipo de elongación la regla de transformación es la 
siguiente. 


la relación inversa es directa 


de aquí que se tenga 





desde luego, las constantes a yb se evalúan a partir 
el nuevo rango de (re. 6.12) valores que se desee 
'ocupenlos pixels dela imagen realizada. Es claro que 
los valores de los pixels de la imagen resizada se 
úObtienen directamente a partir del nuevo histograma 
y de la imagen original 


Otra manera de ver esta elongación lineal y proba: 
blemente más adecuada para efectos computaciona- 
les es la siguiente: supongamos que la imagen origk- 
nalocupa el rango de valores |») y que deseamos 
que la imagen realzada ocupe todo el rango permi 
do |+,.»1) dondels.»! « |; .+yJentonces as fácil vor 
quepormedio de unarelación lineal como!a (6.18) se 
tene 








1) 






histograsa, 
original 


0) 








y 








Fig. 6.3.- El histograma original se transforma de acuerdo a la regla T, porlas reglas de conservación 
discutidas las áreas bajo las curvas se conservan. 
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Fig. 6.4. Ejemplos de elongaciones logarítmica y exponencial. 
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Lo que hace esta simple relación ineal es olongar 
y trasladar la escala de niveles de los pixels para 
Ocupar el rango [vs:1£) 


¡Cuando se desea realzar el contraste de patrones 
Oscuros o claros en una imagen, es conveniente 
"emplear elongaciones logarimicas y axponenciales 
respectivamente, Estas transformaciones, mostradas 
esquemáticamente en la fig. 64, so expresan 
matemáticamente como sigue 


1 + blog(ao) + e 


(6.20 


ca 0 exponencial puede ser necesario un reescala- 
miento lineal (ec. 6.22) para ajustar los valores de la 
imagen realzada al rango permitido. (»;.+<) . 





Otras elongaciones al histograma que son emplea- 
das a menudo son las que siguen. Elongación pro- 
porcional ala frecuencia de ocurrencia de los valores 
que componen los pixels de la imagen 


Voydor1 





(6.25) 








(627) 


grama de la imagen original, a veces se emplean 
Girectamente los valores de sus pixels; esto según 
sea más conveniente para simplificar la transforma- 
ción. Las constantes a, b yc, aligual que en los casos 
anteriores, son para una cierta brillantez 
deseadadelaimagenresultante. La calidad del realce 
se evalúa siempre en términos de una aplicación 
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específica y se hace por medio de la apreciación 
visual de la imagen, desplegada ésta en un monitor 
que generalmente es de alta resolución. En las trans- 
formacionesarriba discutidas, casisiempre serequie- 
re de un reescalamiento lineal. 

La modificación al histograma no necesariamente 
tiene que ser directa, sino que puede ser el resultado 
de una operación o corrección a los valores de los 
pixels d la imagen. Existen situaciones en la investi 
gación experimental en las que el sistema de captura 
de imágenes presenta inhomogensidades espacia- 
les en cuanto a su respuesta ala radiación incidente, 
o bien puede suceder que el haz incidente presente 
variaciones de flujo en su distribución espacial; vea- 
mos algunos ejemplos. En la creación de una neutro- 
grafía ounaradiografía, el haz de radiación empleado 
¡generalmente tiene simerría esférica, por lo que si se 
toma una imagen sin objeto alguno el resultado será 
una distribución que no es pareja, teniendo valores 
más aos alrededor del je del haz y más bajos a 
medida que se incrementa la distancia a ósto, En las 
imágenes de satólte también existo una inhomoge- 
neidad provocada por la atmósfera y la distancia 
variable del sensor alos diferentes objetos (fa. 6.16) 
dela escena. Sin embargo, en este caso es imposible 
obtener una imagen sin objetos para evaluar la inho- 
mogeneidad correspondiente, En cámaras fotográf- 
cas y en vidicones de cámaras co tolovisión, sotienen 
también respuestas diferentes en diferentes puntos, 
debido esto a las aberraciones Ópticas propias de 
estos sistemas, En los casos en donde es factible 
obtener imágenes con y sin escena, se pueden 
emprender correcciones para ajustar la imagen por 
inhomogeneidades del haz, por respuesta variable 

interferencia al sistema de 











captura; 

cia una modificación al histograma de la imagen 
original. Aunque si bien en al caso de correcciones 
complejas como las atmosféricas las modificaciones 
cubren otros aspectos (Lira, Oliver, 1983) de la ima: 
on, Para obtener una imagen que parezca como sl 
hubiese sido obtenida en condiciones de unitormi- 
dad, podemos establecer el siguiente modelo de 
corrección 


806) E, 


6.28) 


que es válido en primera aproximación, En este 
modelo, f(x,y) es la imagon ideal que hubiera resulta- 
dosilas condiciones de respuesta delsistema sensor 
hublesen sido uniformes, g(x y) esla imagenresultan- 
toenas condiciones reales de operación del sistema 
sensor, | (%.y) representa la no-uniformidad de tal 
sistema. Nos interesa determinari(y) para hacer una 
estimación de f (xy) suponiendo que la calibración y 
en consecuencia la respuesta, no cambia en eltiem- 
po de operación y obtención de las imágenes nece- 
sarias para la corrección. Para obtener la función | 
(y) se captura unaimagen de una escena que tenga 
una tonalidad uniforme bien conocida. Para este 
caso, t (x.y) resulta en un campo escalar constante; 


lamémosle k. Sea 9, (4) la imagen real de este 
campo uniforme, de tal forma que se tiene entonces 
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por lo que si suponemos estabilidad del sistema, o 
bien evaluamos | (4,y) a intervalos regulares de tiem- 
po, podemos corregir la imagen por medio de la 
siguiente relación 
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de manera que (xy es finalmente una estimación de 
laimagon corregida. Al igual queen casos anteriores, 
«es posible qua sea necesario hacer un resscalamien- 
to de la imagen corregida. 


Una modificación al histograma que proporciona 
buenos resultados por su aplicabilidad generales la 
que se refiere a la normalización del histograma. En 
estatransformación sele asigna un poso igualaatodos 
los valores de los pixels de la imagen digital; el 
resultado os un histograma plano cuyas barras tienen 
todas el mismo valor, La Imagen respectiva a esto 
histograma plano tiene en principio una relación de 
contraste igual para todas las tonalidades que la 
¡componen y de esta forma se tieno un realce general 
que so puedo aplicar auna gran cantidad de clases de 
imágenes. Aquellas clasos de imágenes donde no 
debe aplicarse esta transformación os donde se io- 
nen imágenes muy homogéneas, o bien donde uno 
varios patrones dominan en ástas, resultando con 
sto un histograma demasiado polarizado a uno.o a 
varios niveles de pixels. En estos casos 

las imágenes realzadas presentan una tonalidad 
tonue y con poco contraste. En casos especificos y 
extremos es mejor diseñar una transformación ade- 
cuada a las necesidades de realce de cada caso., 
Para desarrollar la normalización al histograma, con 
sideremos la integral de la función de densidad. 
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que es desde luego la función de distribución (Rc. 
5.64) y como tal es una función que se increment 
'monctónicamente, Evaluemos ahora la función de 
densidad, o histograma de la imagen realzada, por 
'medio de la transformación dada por la relación 
(6:31). Para esto empleamos la relación (6.17), en 
donde tonamos por substitución directa 
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ya que 





debido alatranstormación que fjamos por medio de 
la integral de la función de densidad (ec. 631). La 
relación (6:32) nos dice que el histograma de la 
imagen realzada tiene un histograma con una densi- 
dad uniforme, que en este caso está normalizada a 
uno. La normalización al histograma es una transfor- 
mación, como quedó demostrado ariba, que no 
depende dela regia inversa T-!,locual es una ventaja 
pues implica un ahorro, ya sea que T-' se calculase 
anallica o numéricamente. Debido a la estructura de 
la transformación del histograma dada por la ecua- 
ción (6.17), era fácil ver que una regia T obtenida a 
partir dela integral de la función de densidad produ- 
ciría un histograma uniforme. Es interesante hacer 
notar esto ya que primero concebimos la normaliza- 
ción al histograma como un realcs que da buenos 





transformación que nos diera como derivada do/dr., 
Precisamente elinverso del histograma original. Ena 
operación computacional de transformación al histo- 
grama, es conveniente tener la normalización depen- 
diendo explícitamente de los parámetros que definen 
a aquél, es decir el número total N de pixels de la 
imagen y el número de niveles L que puede tener el 
rango dinámico de valores de los pixels, en otras 
palabras, L = V,-v, + 1. Entonces en el histograma 
normalizado cada una de las barras tendrá N/L pi- 
xels, si suponemos que éstas tienen un ancho unifor- 
me, de aquí que 


y el ancho de cada barra estará dado por 


por lo que tenemos que el valor de la función de 
densidad relativa al histograma normalizado es 


(626) 


Ahorabien, dela relación (5.15) y, recordando que 
r= (0), tenemos 





por lo que de esta úñima relación obtenemos final- 
mente la regla de transformación para normalizar el 
histograma de la imagen original. 


era 





(6.38) 


Elhecho de contar con una relación que contenga 
explícitamente los parámetros L y N del histograma, 
nos permite evaluar el de salida con los valores reales 
que debe tener en la práctica. En otras palabras, (L - 
1)/N es un factor de escala que ajusta los niveles del 
histograma original para que el de salida tenga la 
densidad uniforme requerida: N/L. En casos reales 
de imágenes digitales es raro encontrar histogramas 
normalizados. planos; debido a la 
naturaleza discreta de los niveles de cuantización de 
los pixels veremos siempre pequeñas variaciones 
que se acentúan sia imagen es pequeña, es decir si 
Ny Lno son números grandes. La relación (6.38) nos 
está dando las nuevas posiciones que van a ocupar 
las barras del histograma de la imagen de salida. 


Las barras del histograma original son desplaza- 
das a nuevos lugares, dentro del rango dinámico 
permitido, en donde dos o más barras pueden so- 
breimponerse. De este manera, reacomodando las 
barras originales, se logra que éstas formen un histo- 
grama aproximadamente plano. Entre más grande 
sea la imagen y más amplio el rango dinámico, más 
fácil será que el histograma de salida tenga una 
apariencia uniforme, Hay que hacer notar que nunca 
sefracciona una barra sino que, entodocaso, las más 
pequeñas se agrupan para alcanzar, lo más cercana 
mente posible, la altura promedio del histograma, 
Como el número de barras es finto y la altura que 
pueden tener es discreta, un histograma perfecta- 
mente uniforme nunca se logra. Otra manera (fg.6:5) 
de verla normalización del histograma, más adecua: 
da para efectos computacionales, es la siguiente: 1.- 
se calcula el valor promedio N/L que deben tenerlas 
barras de salda; 2.- comenzando con las barras de la 
izquierda, se agrupanlas necesarias hasta alcanzar el 
valor más cercano al promedio; 3- las barras así 
agrupadas se asignan al primernivel de cuantización 
de la imágen realzada; 4.- se repite el proceso hast 
agotar las barras del histograma de la imagen origi- 
nal. Es claro que en esta normalización el número de 
barras del histograma de salida es menor al número 
úiginal de barras. En otro algoritmo un poco más 
complejo (Pratt, 1978; Rosenfeld y Kak, 1982) es 
posible alcanzar un histograma plano o casi plano; la 
¡dea aquíes subdividilas barras originales de manera 
que todas alcancen la altura promedio. 


En este método no sólo hay agregación de barras 
vecinas sino también fraccionamiento de las mismas 
para redistribuirsus pixels a vecinos cercanos. Aquíel 
histograma realzado tiene el mismo número de 
barras que el original. Este último método puede 
resumirse en la siguiente relación 
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Fig. 6.5.- Diagrama de la normalización del 
histograma de una imagen con 16 niveles de 





donde laimagen original está cuantizada a K niveles 
y el nivel, tiene pi pixels yen la imagen de salida el 
fivel z, ene q, pixels, de tal forma que 





¡Con MN el número total de pixels y q;= MN/K. El 
"número posible de nivelesz, de laimagen original, es 
igual en principio al de z, de la imagen realzada. En 
algunos casos será necesario elongar el histograma 
"normalizado para que ocupetodo elrango permitido. 
La relación (6:39) debe entenderse entonces de la 

manera: si la igualdad de la izquierda se 
cumple entonces simplemente las primeras k, - 1 
barras son asignadas al primer nivel, 1 de la imagen 
reaizada. Si se tiene la desi Í entonces suf- 
cientes pixels de la barra k, son transferidos a las 
pfimeras k, - 1 barras hasta tener la igualdad, los 
pixels sobrantesson retenidos para la siguiente barra 
dela imagen realzada. Es necesario hacer notar que 
h debe ser suficiente como para que la desigualdad 
dela derecha se tenga al mismo tiempo que la de la 
izquierda. Es posible que la relación (6.39) se tenga 
para dos o más valores de k,. en cuyo caso es 
necesario divicir la barra original para asignar dos o 
más niveles (barras) de la imagen realzada. En resu- 
men, cuando las barras originales son menores a 
MN/Kse agregan incluso en fracciones, cuando son 
mayores a MN/K se subdividen también en trac- 
ciones si esto es necesaro. 

Hasta ahora las modificaciones al histograma se 
han planteado como preconcebidas, es decir, se 
supone a-prior, el tipo de transformación T que de- 
seamos aplicar al histograma original. Esto, si bien es 
conveniente en casos específicos, rigidiza un poco la 
tarea de modificar el histograma; es adecuado por 
anto desarrollar una expresión que permita la inclu- 
sión de cualquiera translormación, desde luego 
dentro de las limitaciones planteadas hasta ahora. A 
esta modificación de tipo general se la conoce como 
especificación del histograma y dados los avances 
tecnológicos puede llevarse a cabo en forma interac- 
tiva por medio de un cursor y una pantalla de desplie- 
gue. Es la generalidad de la transformación la que 
precisamente permite la especificación interactiva, ya 
que la forma del histograma deseado es delineada 
por el cursor en la pantalla, de tal forma que el 
“hardware” del sistema de procesamiento de imáge- 
nes captura los valores que se refieren al histograma 
especificado, estos valores son entonces alimenta- 
dos al "sofware" respectivo el que evalúa latrans! 
mación T necesaria. El resultado de la especificación 
se despliega en la pantalla y se analiza visualmente; 
deser necesario puede repetirse el proceso de espe- 
«ficación. Para expresar matemáticamente la especi- 
ficación del histograma, suponemos que el histogra- 
ima original y especificado se pueden representar a 
través de las funciones de densidad po (0) y Pr(d) 
respectivamente. De aquí que apanir de estas funcio- 
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es podríamos producir sendos histogramas norma- 
lizados, uno a partir del histograma original 


a 


y el otro a partir del histograma especificado. 


ua. [ pa (ula 


(0.42) 


Claramente, las funciones de densidad pu (u) y p 
(uv) son idénticas pues representan ambos 
histogramas normalizados a densidad uniforme, ya 
que el resultado de (6.44) y (6:42) es independiente 
del integrando. Desde luego, larelación inversa « =+- 
6! (udaría los niveles buscados, aunque esto no es 
posibla puesto que la imagen buscada no se conoce. 
Sin embargo recordemos que 7.(u) = puív) detal 
forma que en lugar de emplear los niveles uniformes 
"empleamos los u, por lo que la relación 1=G-! (u) 
proporcionará los niveles correspondientes a la 
función de densidad especificada. Es claro que es 
necesario suponer que G-' existe y puede ser 
obtenida sin mayor ciicutad El procedimiento de 
especificación del histograma puede resumirse 
como sigue: a). normalizar el histograma de la 
imagen original de acuerdo a lo arriba explicado; 
especificar la función de densidad (histograma) 
deseada y obtenerla función de transformación G(9; 
£). obtenerla transformación inversa y aplicar esta 
función a los niveles obtenidos en el paso (a), de 
manera de evaluar los niveles de salida r=G-'(u), de 
esta forma la imagen de salida estará caracterizada 
por la función de densidad Pr(f). En realidad el 
Procedimiento de especificación del histograma 
Puede resumirse en una sola relación matemática, ya 
que a partir de la relación inversa G-! se tiene 





pero los niveles u se obtienen por medio de la 
transformación T, por lo que combinando este hecho 
¿con la relación (6.43), se tiene finalmente 


res] 


Podemos decir que ásta es la expresión para la 
transformación general del histograma a partir de los 
niveles originales. De aquíes claro que la normaliza- 
ción del histograma es un caso particular de laexpre- 
sión (6.44) cuando latranstormación inversa G-! esla 
identidad 6'(1(«)] = 1(c).Al final de cuentas no se 


requiere normalizar explictamente el. histograma 
úriginal, sino más bien evaluar la transformación T 
que es el histograma acumulativo (función de 
distribución) y la función G-! que es el inverso del 
histograma acumulativo de la especificación 
respectiva. Lo que sí es indispensable es poder 
obtenerG-!yasea analticamente en forma detabla. 
Es necesario puntualizar que si bien la formulación 
que se ha hecho para las modificaciones al 
histograma ha sido para variables continuas, la 
versión digitales muy fácil de obtener, Por razones de 


Una de las aplicaciones más importantes de la 
"especificación del histograma es la que se reñero ala 
composición de imágenes para formar un mosaico. 
Alunir dos imágenes paracubrir unacierta escena, 
posible que aparezca una frontera entre ellas dobido 
a diferencias de tonalidad y contraste que pueden 
provenir del hecho de haber sido obtenidas bajo 
condiciones diferentes de iradiación, En ol caso do 
radiografías eltiempo de exposición pudo haber sido 
ilerento, o en el caso de imágenes de satélte que 
hayan sido tomadas en dos épocas del año, las 
condiciones de iluminación no serán as mismas, Con 
el objeto de reducir las diferencias de tonalidad y 
svitarla lormación de una frontera entre las imágenes 
consideradas para el mosaico, es posible modificarel 

de una de ellas para que se aproxime al 
de la otra. Es decir, el histograma de una de las 
imágenes se toma como modelo de especificación y 
a parte de aquí se sigue ol método descrito anterior. 
mente. Otra manera de proceder es definir un tercer 
histograma, arbitrario pero conveniento, tomarlo con 
el especificado y aproximar el histograma de cada 
una de las imágenes a este tercer modelo. Este 
métodoresultatan efectivo que en muchas ocasiones. 
prácticamente desaparecella frontera ontrolasimágo- 
nes del mosaico. En este caso de mosaicos, la espo- 
cificación del histograma recibo el nombre de regis- 
tración de histogramas. En imágenes de satólt, 
puede suceder que hayan variado las condiciones de 
iluminación y las atmosféricas de dos imágenes 
contiguas que se deseen unir para formar un mo- 
saico. En esta situación la correlación de sus respec- 
tivos histogramas no es suficiente para eliminar la 
trontera que se forma en la unión; es necesario 
proceder a una normalización por efectos atmostéri- 
cos, lo que requiere de un modelo de interacción 
atmósfera-radiación. 





8.4.- Operadores Diferenciales 


¡Las modificaciones al histograma planteadas en la 
sección anterior son parte del tema relacionado a 
realces 'de una imagen, el resultado es una 
variación al contraste de la misma pues, de una 
manera u otra, se modifica la distribución de valores 
de los pixels que la forman. 

Otra manera de incrementar el contraste en una 
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imagen es a través de un operador diferencial bi 
dimensional. La razón de esto puede verse al consh- 
dorar la frecuencia de variación de los valores de los 
pixels a medida que nos desplazamos espacialmente 
sobre la imagen digital correspondiente. Si al despla- 
zaros una cierta distancia, encontramos una rápida 
variación, asociamos entonces una frecuencia espa- 
cial alta, de lo contrario la frecuencia que asociamos 
es baja. En forma visual esto se traduce en que 
frecuencias altas estarán asociadas con zonas de la 
imagen donde la tonalidad de la misma cambia rápi- 
damento y las frecuencias bajas corresponderán a 


¡color o gris. Podemos plantear ahora, en el contexto 
do series do Fourier, que una determinada frecuencia 
espacial ostará determinada matemáticamente por 
una función periódica como Senínx) o Cos(nx). En- 
tonces dado que la derivada de Sen(nx) es nCos(nx), 
vemos que entre más alta sea la frecuencia espacial 
(0 más grande) mayor será la amplitud de su deriva: 
da. En forma inversa la integral de Cos(ru) es Sen(n) 
/n por lo que su amplitud será menor para frecuen- 
cias altas. En otras palabras, la derivada de una 
imagon digtal genera un realce de las frecuencias 
espaciales contenidas en la misma, mayor el realce 
entro más ata soa la frecuencia. La integral de una 
imagen, por otro lado, debilita el contraste al atenuar 
la amplitud de las variaciones espaciales; mayor la 

nuación entre más rápida sea esta variación. Un 
¡cambio brusco de tonalidad como el que representa 
sl borde de una clase espectral, un Iineamiento o la 
tronteradeun patrón espacial, estarán caracterizados 
por frecuencias altas. De aquí que un operador dite- 
roncialrosaltará on generallaimagen peroloharácon 
mayor énfasis conos bordes y las líneas dela misma. 
Podemos decir que el operador diferencial es de 
roalco general, pero con una cierta tendencia bien 
definida. Hay que hacer notar que al ruido está carac- 
orizado por frecuencias altas, por lo que si se aplica 
un operador diferencial auna imagen con alto conte- 
nido de ruido el resultado será de pobre calidad. A 
veces es necesario fit la imagen antes de aplicarle 
la operación de derivación. 

Debido al ónfasís con que son realzados los 
bordos ylíneas, a estetipo de operación sale conoce 
como realce de bordes o detector de bordes, aunque 
sl bien la derivada es estrictamente hablando un 
roalco goneral. 

En algunas ocasiones, la derivada seguida de una 
segmentación adecuada, permite retener únicamen- 
to aquellos valores que corresponden a los bordes y 
Iíneas previamente realizados, cono que, finalmente, 
se constituye un fitro pasa altas. En cuanto a la 

misma de derivación, sólo considerare- 
'mos las primeras dos derivadas ya que son las que 
tienen una interpretación directa y sencilla, además 
de que derivadas de orden mayor prácticamente 
convierten a la imagen en puro ruido. 

Las primeras dos derivadas asociadas a una ima- 
¿en g se conocen como el Gradiente y el Laplaciano 
y son un vector y un escalar respectivamente; su 
forma es como sigue, para el Gradiente 





donde i y 3 son dos vectores unitarios en las 
direcciones x, y respectivamente y, para el Lapla- 
ciano 





(6.46) 


Las propiedades de estos dos operadores son 
ampliamente conocidas en cálculo diferencial e inte- 
gral. Existen varias maneras de evaluar estas deriva- 
das: por diferencias finitas, aproximando una super- 
cie a los valores de los pixels de una ventana y 
obteniendo analficamente la derivada de esta super- 
cio y, através de la transformada de Fourier, ya que 





donde ¡ es el número complejo (0,1) y G es la 
transformada de Fourier de la imagen de entrada, Las 
expresiones (5.47) pueden generalizarse para la 
derivada n-ósima. 











donde xj es cualquiera de las coordenadas x O y. 
Demostraremos las relaciones (6.47) empleando 
propiedades básicas del cálculo diferencial e integral 
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de donde se obtiene 





y ya que e* 1 + x, tenemos finalmente 





La respectiva derivada con respecto a y se obtiene 
'nformasimilar alo arriba descrito parala derivada en 
x. La generalización de esto a la relación (5.48) es 
evidente, Con estas relaciones es directo escribirlas. 
expresiones de las translormadas de Fourier del 
Gradiente y del Laplaciano; para el primero se tiene 


(PaCe1)1 e Po, a Je lucy) 





y para el Laplaciano 





De las tres maneras mencionadas para evaluar 
stos operadores, probablemente la más utiizada 
soalade diferencias finitas, dela cual veremos varios 
enfoques más adelante, Como ya se mencionó arr- 
ba, los operadores diferenciales realzan con prete- 
roncialos bordes y las líneas de una imagen, los que 
pueden ocurrir desde luego entodas direcciones, de 
aquí que sea interesante discutir tanto las propieda- 

:s de invarianza frente a rotaciones como las de 
realce direccional del Gradiente y del Laplaciano. 
Como es bien sabido del cálculo diferencial, todo 
operador diferencial parcial: 30/3x*ay9-* — “esun 
operador lineal y cualquiera combinación de éstos es 
también un operador lineal, como los son claramente 
el Gradiente y el Laplaciano. La primera y la segunda 
derivadas representan una medida de la pendiente y 
de la razón de cambio de ésta respectivamente, de 
aquí que podamos considerar a estos operadores 
como locales y por tanto invariantes frente a 
translaciones. Además de estos aspectos generales 
sabemos que el Gradiente y el Laplaciano enfatizan 
las líneas y los bordes los que pueden ocurrir en 
cualquiera dirección en laimagen digital; es deinterés 
por tanto el determinar las combinaciones de 
derivación parcial que generan operadores 
diferenciales invariantes frente a rotaciones. Esto no 
obsta para que en casos específicos sea necesario 
realzar bordes O líneas en direcciones pre- 
establecidas, 





(6.51) 


donde (x.y) yx y” representan los sistemas original y 
rotado respectivamente. Para visualizar las propleda: 
des de invarianza de los operadores diferenciales 
tomemos las primeras derivadas parciales de una 
imagen g(xy) con respecto al sistema rotado (x' 
De aquí que empleando la regla de la cadena dí 
Cálculo diferencial, se tiene 
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(6:52) 
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(0.59) 


por la que claramente se ve que las primeras doriva: 
das parciales de g no son invariantes frente a rota- 
dones. Pero debido a la presencia de las funciones 
Sen0 y Coso en (6.52) y (6.53) es fácil vor que la suma 
de los cuadrados de las derivadas parcialos de y sí 
'ienen propiedades de invarianza. 


(6.50) 


Ahora bien, esta suma de cuadrados es precisa: 
mente igual a la magnitud del Gradiente de g. Enton- 
ceslacantidad |%s/_ representaunaimagen isotró- 
picamenterealzada donda los bordes ylas líneas han 
sido puestos en evidencia con mayor énfasis. 
Veamos ahora que sucede con las derivadas de 
segundo orden; empleando la misma regía de la 
cadena, se tiene que 


Ey AY GA 
Pedro De di e 
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Alsumar miembro miembrolas ecuaciones (6.55) 
y (6:56) se obtiene la siguiente relación de invarianza 


E 





lo que quiere decir que el valor del Laplaciano es el 
mismo en el sistema original que en el rota. Aun 
cuando no lo demostraremos formalmente, ya que 
está fuera del contexto de esta obra podemos, con 
base en las relaciones (6.54) y (6.57), hacer las si- 
guientes inferencias: 1). un operador diferencial lr 
nealisctrópico puede involucrar solamente derivadas 
do orden par y, 2). enun operador diferencial, sctró- 
pico y arbitrario, las derivadas de orden no sólo 
Pueden ocurrir elevadas a potencias pares. En la 
práctica y debido a su fácil interpretación geométrica, 
únicamente se emplean las derivadas de primer y 
segundo orden, combinadas en el Gradiente y el 
Laplaciano. Como ya se había mencionado arriba, se 
pueden diseñar operadores diferenciales que realcen 
nada más en direcciones preferenciales, este tipo de 
"operadores serán portanto no-isotrópicos. Enelcaso 
del Gradiente, al calcular su magnitud y siendo un 
vector, estamos evaluando implícitamente la diec- 
ción de máximo cambio; esto es fácil de visualizar si 
dorivamos a 29/30 conrespectoa +; apani 
de la ecuación(5.51) e igualamos a cero. 


de aquís fácil despejar para ¿con lo que se obtiene 
sl ángulo que forman la dirección de máximo cambio 
con el ejex 





Deestaúltima relación se pueden obtenercos". y 
Sent, para substituir su valor en la ecuación (6.51). 
de tal forma que de relaciones trigonométricas 
simples se tien(6.50) 








(0.61) 


ay 


y que desde luego es la magnitud del Gradiente de. 
La isotropía de ¡"al — se puede interpretar de la 
siguiente manera: toda vez que se presente un borde 
'ounalínea, estaremos calculando su máximo cambio 
independientemente de la dirección que tengan. Por 
oro lado e Laplaciano es un escala y lempro pro 
porcionará el máximo cambi 

Despula de haba acudo ls Prpiadacos de 
invarianzado (%s| yde 1,  ydohaberplantea- 
do su evaluación respectiva por medio de las 
relaciones (649) y (6:50) veamos ahora las 
aproximaciones que existen para dichos operadores 
"empleando el método de diferencias finitas, Una de 
las primeras aproximaciones que se emplearon es la 
substracción del valor de un pixel dado con respecto 
a sus vecinos próximos en las direcciones x, y, par 
esto se definen los siguientes operador 
intermedios 











(0.63) 


(6:64) 


apartir de aquíla magnitud del Gradiente de g puede 
rápidamente obtenerse 





Las dilerencias arriba introducidas pueden ser uti- 
lizadas para estimar el cambio en una cierta dirección 


BeR(103) = Dal 3ICos8 + 8yn(1.3)Seno 


(6.66) 
Esimportante hacer notar que estas diferencias no. 
proporcionan enresicadaluelerde. ll Cen 
el pixel cuya posición es (¡;) sino más bien en la 
posición (+1/271/2), porto que un pequeño error se 
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introduce aquí. Otras aproximaciones centradas en 

(+1/81+1/2),00n + 1/244/2) pueden emplearsa 
Las diferencias de segundo orden 

lenin cier con Secuencia de 

los operadores *, » 3, definidos en (6.63) y (5.54), 

de tal forma que se tiene 


Res) e Cida)) — 281 





y sn forma similar para la segunda derivada en la 
dirección y 





A parti de estas dos útimas relaciones se obtiene 
una expresión aproximada para el Laplaciano de 9 


Las expresiones aproximadas para el Gradiente y 
al Laplaciano de g, están basadas en diferencias 
finitas simples a introducen, desde luego, un error en 
el valor real de estos operadores. Estas aproximar 
clones presentan varios problemas en el tratamiento 
práctico de imágenes digitales. En primer lugar las 
¿iforencias finitas son muy sensibles al ruido, el cual 
está caracterizado generalmente por frecuencias al- 
tas; debido a esto, las imágenes con regular conten- 
do de ruido serán de hecho degradadas si se les 
aplican los operadores mencionados. En segundo 
lugar, los operadores diferenciales así aproximados, 
no son estrictamente isotrópicos, sino que líneas 
diagonales, puntos, esquinas y cambios abruptos de 
dirección de un borde son más realzados que elresto 
de los bordes y líneas que ocurren en una imagen 
digital. En tercer lugar y muy, relacionado con lo 
anterior, las aproximaciones digitales introducidas no 
constituyen tampoco operadores. reversibles, en 
sentido estricto desde luego. 

Otra aproximación al Gradiente y al Laplaciano, 
'queadolece delas res limitaciones quese acaban de 
mencionar, aunque en menor grado, es la que se 
Obtiene al ajustar una superficie a los valores de los 
Pixels do una cia vecindad. Aparte de la expresión 
Aanallicapara dicha superfci, se evalúa (sl y 7 
Para la posición de un porel predeterminado. De 
hecho, para obtener un resultado razonable, es 
necesario ajustar un polinomio p(x,y) de grado malos 
valores de los pixels en una vecindad (ventana) den 
xn, con m = 2 yn = 3 poro menos, aunque m = 1 
yn = 2también se emplean para estimar sl Gradiente. 
Entodocaso elnúmero de coeficientes del polinomio 


¡debe ser menor an. . De esta forma podemos decir 

que el Gradiente o el Laplaciano de dicha superficie 
es la estimación de los mismos en el centro de la 
ventana. Debe precisarse que si bien al aumentar el 
grado del polinomio yal emplear una ventana mayor 
es factible obtener mejores resultados, el cálculo de 
los coeficientes del polinomio involucrado se hace 
más complejo y elaborado. De esta manera y om- 
pleando técnicas bien conocidas del cálculo diferen- 
cial, es fácil vor que el mejor ajuste del plano z = ax + 








Eagrotia) + afaeti Je) = (dd) 41.44) 








lA, del) + (dele) = 620) = 2141.19) 





apartirdelos cuales se puede obtenerla magnitud del 
Gradiente, ya que 


ur 





(6.72) 


El Laplaciano de g no puede obtenerse por medio. 
del ajuste de un plano puesto que las segundas 
derivadas darían un resultado nulo. Si se ajusta u 
supera con un polinomio de segundo orden: x 
+by? + cx + dy + e,con una ventana de 3x 3, el 
desarrollo algebraico es considerable, pero finalmen- 
10 es posible demostrar que el Gradiente tiene com- 
ponentes proporcionales a los siguientes valores 








IS 
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(6.13) 


+31. Jel) y (do. Je1)) 
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(64) 
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'Conelpolinomio de segundo ordenarriba anotado 
síes posible escribir una expresión para el Laplacia- 
no, por lo que después de algunas manipulaciones 
algebraicas se puede demostrar lo siguiente 





Las componentes dadas porlas expresiones (6.73) 
y (6.74) forman lo que se conoce como el Gradiente 
trewit, De las últimas relaciones para y 
"4 podemos apreciar, que al considerar una 
ventana de 3 x 3 cuyo centro es el pixel con 
coordenadas (1), entonces las componentes del 
Gradiente son proporcionales alas diferencias delos 
pxels que se encuentran enlos renglones y columnas 
extremos de la ventana. También vemos que en 
forma similar a la expresión (6.69), el Laplaciano es 
proporcional a la suma de los pixels vecinos al del 
centro menos un múltiplo de éste. Es conveniente 
agregar también que en la práctica no es necesario ni 
adecuado tal ez, el considerar polinomios de orden 
mayor y Ventanas más grandes, ya que podrían 
introducirso inestabilidades en el valor del ajuste al 
polinomio que incrementarían el ruido de la imagen. 





Existon otras aproximaciones para operadores 
diforoncialos de primero y segundo orden que son 
derivadas a part de consideraciones más bien intuk- 
tivas y de emplear otras métricas diferentes a las 
Euclidianas. Aligual que en los desarrollos anteriores 
consideramos una ventana de 3 x 3cuyo pixel central 
es el quetiene coordenadas (¡). Una alternativa alos 
operadores introducidos en las relaciones (6.63) y 
(6.4) es la que se conoce como al Gradiente de 
Roberts, en ásto, on lugar de medir diferencias a lo 
largo delas direcciones x, y, se miden alo largo de las 
agora, poro ques leal lgulert eros. 


eh MleldoJ) = attólJelo 


Intro) — aten 


Otras relaciones que no son precisamente aproxk- 
imaciones a la magnitud Euciidiana del Gradiente, 
pero que tienen una estructura similar, ya que repre- 
úsentan una medida del cambio radiométrico en una 
Imagen digital, son las que siguen y alas cuales se les 
denomina pseudogradiente 


Valp + [864939 — ALI o Its) = 





De estos tres últimos, el que mejor parece respon- 
derala detección de bordes y líneas es el dado porla 
expresión (6.77). Cabe mencionar que estas dos 
últimas relaciones no son deltodoinvariantes frente a 
rotaciones. Es importante puntualizar que lo que es 
lineal isotrópico es el operador diferencial que 
forma de acuerdo a las dos reglas ariba planteadas 
en relación a las expresiones (6:55) y (6:58), 
desde luego no demerita el valor que como detecto- 
res de bordes tienen las aproximaciones aquí discu- 
tidas. El operador de Robers lo escriben algunos 
“autores como una relación de valores absolutos de 








¡Como ya se había mencionado, los operado: 
diferenciales resaltan también ol ruido, sobre todo, 
ésto es de alta frecuencia, de tal forma que silarazón 
señal/ruido no es adecuada entonces la imagen 
transformada puede resultar firalmente degradada. 
De aquí que sea conveniente en ciertos casos aplicar 
operadores formados más bien por diferencias de 
promedios; de entre éstos el más conocido es proba- 
blomente el de Sobel 








Sinembargo, estetipo de operadores diferenciales 
como el de Sobel producen un ensanchamiento de 
los bordes resizados lo que para algunas aplicacio- 
nes como el seguimiento de líneas resulta inadecua- 
do. Por otro lado, debido a los promedios introduci- 
dos, se realiza una operación de firaje, que suaviza 
el ruido pero que hace en principio ireversible la 
"operación de diferenciación. Existen otros operado- 
res para detectar bordes y Ííneas pero no son de 
naturaleza diferencial, aunque al igual que éstos 
pueden ser establecidos como una operación de 
convolución entre una ventana con factores de peso 
y la imagen de entrada. También hay variantes de 
estos detectores en donde se introducen relaciones 
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no+lineales; este tipo de operadores serán revisados 
en el capítulo sobre fitrajes. 


6.5. Transformada de Fourier 


Uno de los exponentes más importantes de las 
transformaciones que pueden hacerso a una imagen 
os la de Fourier; sus propiedados y aplicaciones han 
resultado muy interesantes en el estudio de estructu- 
ras bidimensionales y en el desarrollo de fitrajes 
Gigtles Ora propieads importantes serán tam- 
bién revisadas en esta sección, proporcionando una 
visión general de esta transformación, ya que proba- 
biemente sele ha dado un énfasis demasiado grande 
aldisoño de fitros por medio del espectro de frecuen: 
clas que genera la transtormada de Fourier. 
Latransformada de Fourier es una transformación 
lineal, invenible y con kernel separable, por lo que 
Puede serutlizadala relación (5.6) para introducir una 
dofnición adecuada. De esta forma el komel de la 
transformación de ida es el siguiente 








para una imagen cuadrada de dimensión N. Elkerel 

de la transformación de regreso es simplemente el 

conjugado complejo de. Es fácil comprobar que los 

kornols le ida y de regreso de la transformada de 

Fouror satisfacen las relaciones (6.7). Lo de arriba 

Quiere decir justamente que la transiomada de 
'ourler discreta es la transtormación Z 


sobreloscírculos unitarios paralos cuales es válidolo 
siguiente: 


ed espa 





> explain 


y son precisamente las propiedades de simetría de 
stas n-ralces de la unidad las que hacen posible la 
elaboración de algoritmos computacionales de muy 
rápida ejecución. Sinos referimos ala relación (5.2), 
vemos que para la transformada de Fourier, las mat 
ces (:) y lo) sontales que 


LPI = [hoy 7 LO) = [8 


donde en general los elementos (m,n) de la matriz 
están dados por 


1 
hexol2n30m/X). con mn 





5) 

para una imagen cuadrada. La transformada de ro- 

greso está definida por las matrices inversas ||, 
y [5izz  cuyoselementos (4,1) en general 


PEI MLN). com Ral La daa 


y nótese que el factor de normalización 1/N se ha 
Puesto en latransformación deida. Con esta formado 
ver las cosas es claro que 


donde es la matriz identidad, ya que las ex: 
ponenciales involucradas cumplen con la siguiente 
propiedad: 


lopara ko 
ol 21 00/5lexpl 29 3m0/1) > 
park 


0) 
La forma discreta de las exponenciales en la defini 


ción de los kemels A y B implica que estas funcion 
son periódicas 





pl d/31 = aspl20K (040/31 = oxpl27Jo(Kon)/8 


6.40) 


ya que esta función exponencial puede ser escrita 
Somo una combinación compleja de funciones Seno 
y Coseno las que cumplen con propiedades de perio- 
“cidad. Esto quiere decir que la secuencia 0 imagen 
de entrada, ((»,.1.)) - yla secuencia o imagen de 
salida, (7(%,.x,/) _ son secuencias periódicas, ya 
que la transformada discreta de Fourier (TOF) 
define. como 





E Eran dorpl Pr (sm da0)/81 





(6.00) 
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para una imagen cuadrada. Es decir la TDF de una 

imagen es esencialmente una representación en 

series de Fourier de un campo escalar bi-dimensio- 

únal; para que esto sea válido, el campo debe ser 
ico. 


La TOF bidimensional posee una serie de propie- 
dades de simetría muy interesantes para el análisis de 
imágones, las consideraremos sin demostración 
matemática ya que existen tratados especializados 
para esto. Las propiedades de uso probablemente 
más extendido son las siguientes 


Estructura de una Imagen 


Debido a que la TDF de una imagen es finalmente 
representación en series de Fourier de la misma, 

. fácil ontender que esta transformada proporciona 
ladistribución de frecuencias delas variaciones delos 
valores de los pixels que componen la imagen. En 
otras palabras la TDF es una descomposición del 
campo bidimensional en un espacio de frecuencias 
'ompleando como funciones base al Seno y al Cose- 
no. Las frecuencias a las que nos referimos son 
entonces la rapidez a la que cambian los valores de 
los pixels en ambas direcciones del espacio de fre- 
cuencias; de aquí que la TDF permita la evaluación 
del espectro de frecuencias espaciales de una ima- 

Esto ospectro de frecuencias o distribución de 
variaciones delos valores de laimagen, enrealidades 
una medida de su complejidad, La TDF es por tanto 
un indicador de la estructura de una en el 
pS ee sole introductoria de este 
trabajo, 

De osta manera la distribución de frecuencias 
puede ser empleada para dilucidar ciertas propieda- 
des geométricas de una imagen; por ejemplo, las 
Iíneas de drenaje obtenidas a partir de una imagen de 
satélite, tendrán un espectro de frecuencias caracte- 
rístico de la clase a la que pertenecen. En otras 
palabras dos patrones de drenaje mostrarán espec- 
tros similares si pertenecen a la misma clase de 
drenajes, aun cuando sean aparentemente diferentes. 
porvariaciones en suorientación y extensión. De aquí 
que la TDF ponga de manifesto o realce, ciertas 
Propiedades geométricas de periodicidad de los va- 
lores de una imagen que no están evidentes en el 

espacio original de coordenadas. En resumen, imá- 
genes con propiedades globales de periodicidad 
similares, tendrán espectros de frecuencias pareci 
dos. Además de esto, los espectros de frecuencias 
son más fáciles de comparar entre sí, que tratar de 
visualizar similtudes geométricas en las imágenes 
originales. La TDF es portanto una representación de 
laimagen digítal en un espacio donde se encuentran 
roalzadas características espaciales de la misma; por 
esto, la TOF puede pensarse como un realce global al 
campo bidimensional de entrada. 











Teorema de Shiting 


Sila secuencia (Xik.ks)) + Kits La 
la TDF de la secuencia <(1;.0,)). 0300 = 1.2, 
N, entonces la secuencia 





tiene la siguiente TDF 


lo que quiere decir que un cambio de posición en el 
espacio de coordenadas implica un cambio de fase 
en el espacio de frecuencias y viceversa un 
úcorrimiento en frecuencia trae un cambio de fase en 
el dominio de la imagen, entonces la secuencia 
Xky-ks kk) corresponde a la TDF de 


epldrjóhon 4 Kon 1ÚNIC nin 


donde k, es un desplazamiento arbitrario en 
frecuencia. 


Teorema de Convolución Circular 


Aesteteorema se lellama circular por la propiedad 
de periodicidad arriba expuesta y para diferenciarlo 
del teorema líneal en donde no se suponen 

dice lo 





Aida AX) 





¡son respectivamente las TDF de las secuencias 





entonces la TDF obtenida del siguiente producto 


Uh) Me Maeda PK 


corresponde a la secuencia! x,(s,...)Jestá dada por 


(6d 


Este teorema es análogo al correspondiente para 
funciones continuas, sin embargo dada la naturaleza 
periódica de las secuencias 'x:im.n001 y (otm 
2. cuando una secuencia es desplazada más allá 
del último punto de la otra secuencia (al realizarla 
convolución), no encontrará valores cero como 
debiera, sino la repetición periódica de la secuencia 
en cuestión. Entonces para forzar a la convolución 
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discreta a que tenga los mismos valores numéricos 
'uelaconvolución lineal de dos secuencias finitas, es 
necesario construir dos secuencias aríficiales que 
tengan valores cero en los extremos. 


Cantidades Escalares Relacionadas a la TDF 


Latransformada de Fourier de una función real, es. 
¿generalmente compleja, es decir 


(Ko) RR ke) 6 TA 


donde R es la parte real e ' la imaginaria, a partir de 
aquí se definen las siguientes cantidades 


-El 





¡pectro de Fourier o de trecuencias 





Estas tres cantidades son las que 
cuantitativamente nos dan una visión general de la 
estructura de una imagen en el sentido arriba 
expuesto. 

A diferencia del caso continuo, no es necesario 
preocuparse acerca dela existencia dela TOF, ya que 

existe siempre en el caso discreto 
puesto que la imagen de entrada siempre está 
acotada. Además es importante hacer notar que la 
TOF noes una aproximación del caso continuo, sino 
que es una transformación completa. La gran 
mayoríadelos espectros de frecuencias de imágenes 
decaen rápidamente como una función de la 
frecuencia y de aquí que los términos de frecuencia 
altatengan la tendencia a verse oscurecidos cuando 
son desplegados para su inspección visual; es decir, 
el lóbulo central es mucho mayor que los lóbulos de 
segundo, tercer y orden mayor, por lo que en un 
sistema de despliegue que tenga disponible 256 
niveles de grs (o de color), al lóbulo central le 
Corresponderá eltono relativo al 256 y los lóbulos de 
segundo y tercer orden tendrán tonos muy oscuros 
cercanos al cero. Debido a esto se emplea una 
técnica adecuada que compensa los efectos de 
despliegue causados por esta particularidad del 
espectro de Fourier, estatécnica consiste en utiizarla 
función 








enlugarde '//%u.K:11  .Desdeluegolos valores 
de Dik,k,) deben estar normalizados a la escala de 
despliegue disponible. El uso de la función Mk: +) 
preserva los valores cero en el plano de las 
frecuencias ya que DMti.ka) +0 cuandoli( ke 
También hay que hacer notar que 1%.) 
es una función no-negativa y que se emplea. 
únicamente para efectos de despliegue de la TF, 
Cabe también mencionar que para estos efectos y 
para algunos cálculos, los cuadrantes de la TDF son 
reordenados para que el lóbulo de orden cero 
aparezca en el centro de la imagen resultante, Est 
cuadrantes pueden reordenados 
automáticamente multiplicando a la imagen original 
por el factor (+)! antes de ejecutar la TOF, 
Por lo que debido a la identidad 






la TDF puede reescríbirse como. 





Soparabilidad 


En la definición de la TDF ya se anotó que ésta es 
una transformación con kernel soparablo, lo que 
Quiere decir que la principal ventaja de la propiodad 
de separabiidad es que la TOF puede obtenerse por 
dos aplicaciones sucesivas de la transformada uni- 
dimensional, es decir 


Jm 





ñ 
¿a k1) e $ E tna Jexpl-2r)K 203/81 
0) 


Sin embargo, en la práctica, implantar la TDF de 
esta manera implica el uso de muchas operaciones 
de entrada/salida de la unidad de almacenamiento 
donde resida la imagen original, de contarse con 
suficiente memoria central o con un procesador de 
arreglos, es mucho mejor emplear la versión bi- 
dimensional de la TDF. Para hacer más eficiente la 
evaluación de la TDF se toman imágenes con 
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dimensión de orden N = 2? con n entero, con lo que 
se aprovechan las propiedades de simetría de las 
raíces de la unidad, obteniéndose de esta manera lo 
quese conocs como transformada rápida de Fourier. 
Las raíces o valores del keel en cuestión pueden 
estar pro-almacenados en tablas para una serie de 
valores de n y si la imagen no es cuadrada de orden 
2, se puede entonces agrandar con unos, hasta que 
tenga esta dimensión, el resultado no se altera 
básicamente ya que con esta técnica sólo se agrega 
una frecuencia constante a la imagen original. 


Perlodicidad y Conjugado Complejo 


La TDF y su inverso son periódicas con período N: 
Teorema de la Convolución Circular, por lo que se 
tiene. 


AA) 


donde m y n son enteros y N la dimensión de la 
Imagen. Por otra parte 


SN) 


o bien 
[70431 + [Fa 0] 


Dado que la TDF es en general compleja, el 
dominio de frecuencias tiene por tanto 2N 
componentes, incluyendola parterealyla imaginaria, 
por esto tal vez podría pensarse que esta 
transformación produce un incremento en la 
dimensionalidad de la imagen; no es esto el caso sin 
"embargo ya que debido a la propiedad (5.96) la casi 
mitad delos tórminos dela TDF son redundantes, por 
lo que la dimensionalidad se conserva. 





Rotación 


Sigín,n,)esrotada porunángulo 1. ,entonces 
la respectiva TDF es rotada también por el mismo 
ngulo.En forma análoga, larotación de /(x,.+; trae 

consigo la rotación de g(n,.n¿) por 
«Imismoángulo.Sinembargo, es necesanotortaren 
¡cuenta quel rotar unaimagen, los lugares asociados 
alos pixels rotados no necesariamente coincidirán 
¡con los lugares de la malla o latice empleados en la 
imagen original, de tal manera que será necesario 
realizar una interpolación como la convolución 
cúbica. Los valores del pixel original más cercano al 
rotado, así como los próximos vecinos de aquél 
contribuirán al valor correspondiente del pixel de la 
rotada. Además es necesario tomar en 

cuenta que la imagen rotada y la original serán así 





siempre referidas a la misma lattice. 
Distributividad y Escalamiento 


Dado que la transformada de Fourier es lineal, se 
tiene 


O 


y on general 


Fis (ass Julian] A 7140 (a 10031 03) 


(oo 
Además de esto si a y b son escalares 


Estas relaciones se siguen claramente de la 
definición misma delatransformada de Fourier, Dado 
que las variables n,, n,. kk, y los valores que puede 
asumirla imagen son discretos, es necesar 
una interpolación en forma similar al caso 
rotación. 






Valor Promedio o Componente de Frecuencia 
Coro 


Si tenemos el siguiente valor promedio de una 
imagen 


O 


y si adomás consideramos la componente de 
frecuencia cero de la transformada de Fourier de la 
misma imagen 
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entonces se tiene que 
2-70.) 


(6104) (6.100) 
Enanalogíaconlateoría electromagnética 710,0) 
sele conoce también como la componente d.c.de la 
imagen yen muchas ocasiones es restada al espectro 
de Fourier para faciltar su análisis, 





Gradiente y Laplaciano 


Yase anotaron arriba las relaciones para obtener el 
Gradiente y el Laplaciano por medio de la 
transformada de Fourier, también se incluyeron las 
rospectivas demostraciones. Estas relaciones son 
muy importantes para evaluar dichos operadores 
iforenciales los que son ampliamente utiizacos para 
rosaltar las líneas y los bordes de una imagen. 


Convolución y Correlación 


Para sistemas lineales, invariantes frente al tiempo 
sotione, que si fy g son dos imágenes cualesquiera 
y si F y G son sus respectivas transformadas de 
Fourier, entonces 


y además 


Lor y e ros 


dondet* geslaconvolución de fy g como yase anotó 
arriba y donde 


£loi05) o glam) EJ ta. 8atnrta, 0348) 


EJ 
sé 


(6.105) 


es la correlación entrefy g. La validez delteorema de 
la convolución es la que permite la ejecución rápida 
de fittajes en el dominio de frecuencias y faciita 
también la comparación entre dos imágenes o entre 
dos patrones, así como la estimación de la textura a 
través de la función de autocorrelación. Cabe 
mencionar también que la convolución de una 
imagen con un templete puede constituirse, bajo 
ciertas condiciones, en un fitraje lineal a la imagen. 
Esta técnica puede generalizarse para una clase de 
fitros no-lineales y detectores de líneas y 
ramificaciones. El tema de fitrajes será desarrollado. 
en el próximo capitulo. 

Existen otras transformaciones similares en 
naturaleza a la de Fourier (ver figuras 6.6a- 6.59). Las 
transformadas de Walsh, Hadamard, Coseno, Haar y 


Sant, son también unitarias y desde luego de kernel 
separable. — Constituyen — asimismo — una 
escomposción del imágen cial en funciones 
base (fg. 6.5) de un espacio generalizado de 

frecuencias; son diferentes a la de Fourier en el 
sentido de que fueron concebidas en primera 
Instancia (con excepción de la de Coseno) en forma 
digtal, mientras que la de Fourier fue definida 
Primeramente para señales. También su 


mismo son empleadas para la comprosión de 
imágenes que requieren de ser almacenadas o 


Latransformada Walsh para una imagen digital y ae 
define como sigue 


o 
Dat AA 
o o o 10 





009, 
TN 


(6.106) 


¿onde se supone que eltamaño de la imagen es de 
y q Py donde bm) es elbásimo bit 
delarepresentación binaria de m, por ejemplo sim = 
6 (110 en binario), entonces bm), = m1 y 
Dl =$ dodo luego EN 

necesario calcular pH y después evaluar 
lareprenertación rara Cometrondent, Eltórmino 
1/N es un factor de normalización asociado al kernel 
de la transformación de ida. La transformación 
inversa tione el mismo kemel salvo por el factor de 
normalización que en este caso es 1/M. Nótese que 
silaimagen es cuadrada entonces M = Ny los komels 
de Js rnefomacions de da, y teeo son 
idénticos. Hay que hacer notar también que 
Gerencia dela transformada de Fourier la de Wal 
esreal porlo que requiere menos memoria digital y es 
más rápida. 


La transformada de Hadamard es muy parecida 
(69. 6.5) a la de Walsh. Cuando las imágenes 
procesadas son cuadradas y de orden 2, las 
transformaciones de Hadamard y de Walsh tienen las. 
mismas funciones base aunque en diferente orden. 
Es por esto que en la iteratura, se habla 
frecuentemente de la transformada Walsh- 
Hadiamard. La transtormada de Hadamard se define 
de la siguiente manera 


2 


CN 
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Fig. 6.6.- Funciones base de las transformadas de Fourier (a,b), Hadamard y Walsh 
(e), Coseno (a), Haar (e) y Slant (9. 
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Donde, qyb/(m) son definidas de la misma forma. 
que en la transformada Walsh. Nóteso que se tiene 
tambión un factor de normalización 1/N, En forma 
análoga ala transformada de Fourier, es conveniente 
expresar la de Hadamard de tal manera que los 
argumentos del kernel correspondiente se 
incrementen a la par que la frecuencia de las 
funciones base asociadas (fg. 6.6). Por esta 
razonamiento, la de Hadamard se convierte en una 
transformada através de la cual descomponemos en 
frecuencias a la imagen original. Esto ordenamiento 
particular se logra reescribiendo la transformada de 
Hadamard de la siguiente manera 




















dondes.(u) = b,, (u),5,(4) = b,,,(u) + b, (4) y, 
= bu) + b(U). Con lo Gual vemos q 
translormada ordenada de Hadamard 
ordenada de Walsh para N 












lastransformadas Walsh y Hadamard es claro que el 
mismo algoritmo puede servir para calcular los 
Komnels de iday de regreso, Dado que latranstormada 
ordenada de Hadamard puedo verse como una 
descomposición en frecuencias de la imagen 
original, podemos pensar que proporciona 
información acerca de la estructura de la imagen en 
forma análoga a como lo hace la transformada de 
Fourier; lo mismo se puede ceci de la de Walsh, 
Otratransformacadeuso generalizado, sobretodo 
para la compresión de imágenes, es la transformada 
Coseno, a partir de la cual y bajo ciertas restricciones 
(ergodicidad es posible evaluar la descomposición 








decuentas, una modificación de la versión discreta de 
la de Fourier. Está desarrollada a partir de funciones 
base compuestas por funciones trigonométricas 
'Coseno de fase cero y de naturalezareal y discreta; se 
define como sigue 


E atnon)(Cos(2n01 lr Cos(2041)10] 





El kernel de esta tranformación consiste 
básicamente de funciones Coseno muestreadas. La 


transformada inversa es iguala de ida excepto por 
Un factor constante y uno de normalización 


a cbraLico 
mi 0 





¡donde el factor de normalización se ha colocado aquí 
nlatranstormada inversa y donde el factor constante 





Mal a 





con lo que se ve que el primer término de (5.110) es 
el valor promedio de la imagen de entrada. Las 
transformaciones de Haa y Slant son muy similaros a 
las ya descritas en esta sección por lo que no se 
"estudiarán explícitamento, sin embargo, a man: 
Ilustración se muestran Sus respectivas funciones 
base en la figura 6.6. 








6.5.- Descomposición en Component 
Principal 








Las transformaciones a la imagen que homos visto 
hasta ahora son aplicables a imágenes mono- 
espectrales, o a cada banda por separado de una 
imagen muliespectral, La descomposición en 
Componentes principales es una transformación al 
conjunto de datos que se aplica simultáneamente a 
todas las bandas que componen la image 
salida se obtiene entonces una imagon mul 
también, Esta transformación se conoce pol 
autores como Karhunen-Loeve (KL), Hoteling o 
Transformación de Eigenvector. En el pre 
trabajose prefierela denominación anotada paraesta 
sección debido a la interpretación geométrica que 
eta tlono, ya que la imagen de salida se encuentra 
un espacio generado por vectores base 
cuyas direcciones son precisamente las de las 
¡componentes principales de la imagen original. Esta 
transformación es desde luego unitaria y por tanto 
puede considerarse como un realce global a la 
imagen tal y como se verá en el desarrollo 
subsecuente, Tiene la gran ventaja, que aligual quela 
transformada de Fourier, de no requerir parámetrc 
de entrada, sino que hace uso nada más de las 
propiedades estadísticas del conjunto de la image 
Multiespectral. Además caca banda de salida puede 
toner una interpretación fi que los 
sigenvectores y elgenvalores asociados la tienen en 
principio. Otras propiedades se discutirán através del 
desarrollo matemático que involucra esta 
transformación. 

Veamos ahora una definición formal suponiendo 
¿que la imagen de entrada es un campo aleatorio, 
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resultado de un proceso estocástico generalizado de 
la maneracomo se discutió en el capítulo precedente. 
Después de esta definición haremos una 
reconsideración de la transformada aportando 
elementos físicos y geométricos, Consideremos la 
clase deimágenes formada por y bandas donde cada 
banda! (+) representa una imagen digtal resultado 
de la realización de un proceso estocástico. Sea (1) 








El sontido vectorial de (5.111) consiste en que el 
operador A actúa sobre los valores del campo que 
lorman los pixels de la imagen multiespectral de 
entrada, para producir el campo vectorial g de la 
misma dimensión que el original. La transformación 
(6.111) está sujeta ala restricción de que la matriz de 
covarianza K, de la imagen g es diagonal. En el 
'spacio de shlida la matriz de covarianza es por 
definición 


Larazón de ser de estarestricción sobreK, consiste 
en lo siguiente: En general existe una cérrelación 
ntrolas bandas de una imagen mukiespectral, sobre 
odo cuando éstas son contiguas, ya que una banda 
+s adquirida para un cierto intervalo de energías, es 
de esperarse entonces que para intervalos cercanos 
o próximos la infomación correspondiente sea 
similar. En otras palabras entre bandas contiguas 
habrá una redundancia de información, que se 
traduce en que la matriz de covarianza de la imagen 
úriginaltiene elementos diferentes de cero fuera de la 
diagonal. Imponer una matriz de covarianza diagonal 
para la imagen de salidaimplica que se elimina dicha 
redundancia de información o correlación entre 
bandas, lo que tras consigo una reducción en la 
dimensionalidad de la imagen mutiespectral. La 
matriz A de la relación (6.111) es una matriz untara, 
cuyos renglones están compuestos por los 
vectores de la matriz de covarianza K, de las 
bandas originales. La matriz A realiza una 
diagonalización de la matriz K, de tal manera que la 
matriz do covarianza de la imagen transtormada 








es una matriz diagonal cuyos elementos son los 
sigenvalores de K, arreglados en forma decreciente. 
A la relación (6.113) se le conoce como una 
transformación de similaridad (Arfcen, 1970) yA esun 
operador hermitiano y real. La transformación KL. 


produce un conjunto de bandas decorrelacionadas 
que forman una imagen con dimensionalidad menor 
ala original, Este proceso requiere, desde luego, del 
conocimiento de la matriz de covarianza definida por 
las bandas de laimagen mutiespectral original. Dicha 
matriz debe ser modelada, estimada o medida según 
sea el caso; en este trabajo veremos un ejemplo de 
modelación. En forma explíctala transformada KL se 
escribe como 





aun) + EE IIA AO 


3% 


(6116) 


donde el keel A satisface la siguiente ecuación de 
eigenvalores y elgenvectores 


OA) 


(6115) 


donde » representa los ejgenvalores dela función de 
covarianza y A los sigenvectores de la misma. En 
otras palabras, la matriz detranstormación Asatístaco 
la relación 


MM 


(6.116) 


Es importante hacer notar que la matriz K, es 
simétrica y por tanto diagonalizable (Tyidesloy, 1975) 
con eigenvalores reales. Es decir, que los 
eigenvalores 4 deben cumplir con la ecuación 





secular 
Ig - a] =0 


(000) 


donde! esla matriz identidad. Los eigenvectores 
és — correspondientes la ecuación (6.115) 
satisfacen la relación 


Kg + dada 
Aa 


(6.118) 
estos eigenvectores forman una base ortonormal: 


Cabe agregar que se cumplen las siguientes 
relaciones 


(6.119) 


donde “1, es la dela de Kronecker. A los 
'igenvectores se les llama componentes principales 
y como los eigenvalores están arreglados de manera 
decreciente, entonces el contenido de información de 
la imagen es redistribuido entre estas componentes 
(o direcciones principales), de manera que la mayor 
información está concentrada en las primeras 
componentes o bandas de la imagon de salida, 
mientras que en las restantes se tiene muy poca 
información y prácticamente todo el ruido. 


Veamos ahora, en forma gráfica, cómo podemos 
visualizar a la descomposición en componentes 
principales. Supongamos que tenemos una imagen 
de tres bandas por lo que en un espacio 
tridimensional (fg.5.7) se puede representar el total 
de datos en forma de nube. Esto se logra graficando, 
elvalor quetienecada pixelenlastres bandas, eneste 
espacio tridimensional, el cual está construido por 
tros ajos normales entre síy cuyas respectivas escalas 
sonigualos alos rangos dinámicos que pueden tener 
lospixels encadabanda. Asíun pixel cuyo valor enlas 
tres bandas sea, digamos, (63, 65, 72) tendrá 
justamente las coordenadas (63, 6S, 72) del 
mencionado espacio; estas consideraciones son 
fácilmente generalizadas a imágenes de bandas. El 
conjunto de datos gralicados de esta manera, forma. 
lo que se conoce como una nube de datos y que se 
muestra esquemáticamente en la figura 6.7 (ver 
tambión figura 5.2). Los datos así representados son 
toforidos a un espacio de tantas dimensiones como 
bandastengala imagen, en el presente ejemplo, en al 
jo x están los valores de la banda 1, on ely los de la 
banda2 y enelejezlosdela banda3. Sila correlación 


datos presente una elongación que dencte dicha 
correlación. En la figura 6.7 se puede apreciar que 
entre la banda 1 y la 2 existe una alta redundancia de 
información.Por esto os fácil ver que si los ejes 
Orgnales (o. beb) sonrtados, on inclusión de 
Unatranslación delorigon, detal maneraqueel ejeb,, 
se coloque a lalargo de la máxima elongación de la 
nube de datos, el b, a lo largo de la siguiente 
elongación en da y así sucesivamente, 
tendremos entonces un nuevo sistema de 
coordenadas (b',, b, b'). Sobre el ejo »; queda 
definida la dirección de mayor elongación de los 
alos y que puede ser obtenida, por ejemplo, por 
mínimos cuadrados. En un sistema de tres 


sólo podremos reorientar secuencialmente el eje b, 
yaqueelb, quedará automáticamente reorientado ya 
¡que hemos supuesto un sistema ortonormal; aun así 
elejemplo es suficiente para nuestros propósitos. Las 
ecuaciones que describen está mutirotación de ejes 
son las siguientes. 








Fig.6.7. Nube de datos de una imagen multi 
poctral de tres bandas, el vector de posición 
mostrado coincide para la imagen original y 
para la resultante de la transformación KL. 
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donde *i_ son las direcciones del nuevo espacio y 
representan los valores de los pixels do la imagen de 
salida para cada banda de ésta, b, son los valores de 
entradaya,, son los coeficientes dé latranstormación. 
Esta Úlima relación nos dice que básicamente la 
descomposición en componentes principal 
transformación reversible, que seforma por medio 
una combinación lineal de las bandas originales 
tomando el valor de cada pixel como entrada al 
proceso, Los coeficientes a, se determinan vistas las 
cosas así, porel grado de! aplicado alos 
del espacio original. Por medio de esta 
transformación puedo verse fácimente (fgura 6.7) 
Que elintervalo dinámico de valores para elnuevo ejo 
b, es mucho mayor que el original, que el intervalo 
Para b, es mayorque elb; y que el que se refiere lb, 
es mucho menor que el b,. 

Enotras palabras, la mayor parte de la información 
hasido puesta sobre el nuevo eje b; (banda de salida 
b;),una fracciónde estasobre elb;y casinadao puro 
nudo sobre el b,. En esto ejemplo, es claro que sila 
correlación entre bandas es muy alta, la 
“imensionalidad de la imagen de salida puede 
reducirse a casi uno. Para las imágenes Landsat 1,2 
y 3 de cuatro bandas y cuya correlación puede 
sobrepasar el20%, sotiene que después de aplicarla 
descomposición en componentes principales la 
imensionalidad se reduce a aproximadamente 2.2; 
aestosele conoce como dimensionalidad intrínseca. 
Para los Landsat 4 y 5 de siete bandas, la 
Gimensionalidad intínseca es de entre 3 y 4, Cabe 
anotar aquí que a las nuevas direcciones se les 
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¡conoce como ejes o componentes principales y por 
sto el nombre de esta transformación. 
Hasta ahora hemos introducido la transformación 
KL a través de una definición formal, 
matemáticamente hablando y hemos. discutido 
conceptualmente las consecuencias que trae 
consigo empleando un esquema gráfico (fg. 6.7). 
También algunos aspectos físicos fueron 
mencionados, como es el hecho de la redundancia 
en información de bandas contiguas de la imagen 
'multiespectral, Esto proviene de la consideración de 
que pensamos, que, en general, el comportamiento 
'sperado de un sistema bajo estudio, debe ser muy 
lar cuando se observa en intervalos de energía 


contiguos, 

Para fines computacionales y para profundizar aún 
más en esta transformación, revisemos con mayor 
detalle la definición hasta ahora formulada. 

En primer lugar, es necesario entender qué se 
quiero decir con latranstormación detipo vectorialde 
la relación (6.111. En esta operación la matriz A 
actúa sobre cada pixel de la imagen de entrada, 
tomando aóste como un vector contantos elementos 
como bandas tenga la imagen original. Elvalor de un 
pixel de entrada será entonces t(1), cuya localización 
rolativa en la imagen estará definido por su vector de 
posición r; f es por tanto una función mutivaluada 
paratodas las bandas de la imagen mutespectral. De 
sta manera el operador A transforma cada pixel de 
entrada para producirlos pixels de salida que forman 
asu vezlaimagen resultante y(s), cuya dimensiónes, 
estrictamente hablando, igual a((r), pero con bandas 
decorrelacionadas, porlo quela dimensión intrínseca 
de g(s) es menor a la dimensión física de (0. 
Entendiendo esto veamos con detalle cómo se 
roaliza el cálculo y bajo qué condiciones se evalúa la 
transformación KL. Sabemos que cuando setiene un 
campo aleatorio homogéneo y bajo ciertas 
condiciones estadísticas (ergodicidad) la media o 

Primer momento de la imagen mutiespectral 
Conespendare se puede estimar con la siguiente 
relación 








donde : denota valor esperado, 1, es el vector que 
representa los valores que componen al pixels + 
ésimo de la imagen original a través de las bandas 
que la componen y N es el número total de pixels de 
la imagon; nótoso que en este orden de ideas 
también es un vector. A partir de la relación (5.122) es 
fácil obtener la matriz de covarianza de la imagen 
original. 


e E A 0 


6.123 


dondet denota transpuesta de la matriz columna 

(Es up, ; Bsjolas miamas euposiciones quese 
hacen para la media, es factible encontrar una 
expresión aproximada para la matriz de covarianza 





(6,124) 


Es claro que la matriz de covarianza es de 
dimensiónigual alnúmero de bandas de laimagen de 
entrada y además es una matriz simétrica ya que la 
correlación de la banda 1 con las + 1, es igual a 
¡correlación de la banda k con la 1. La correlación de 
una banda con ella misma es, desde luego, igual 
uno. Cabe anotar que los elementos de la matriz de 
Correlación (£) = (2. in están dados porla siguiente 











(6.125) 


dondelosíndices my nvan de 1 alnúmero de bandas 
delaimagen multiespectral + “=.  sonlos elementos 
dela matriz de covarianza! +). La relación (6.125) 
válida tanto para la imagen de entrada como para la 
de salida. En la práctica no es necesario calcular la 
matriz de correlación para realizar la transformación 
KL; esta matriz se obtiene en este caso nada mas 
¡como una manera de estimar que tanta redundancia 
sn información existe en la imagen original, Ahora 
bien, la media de la imagen de salida está dada por 





Uy + <leto)] EAN Y Ae 1G0)] = A 


(6.126) 
esto es posible ya que [1] es una matriz con 
elementos constantes y definidos para una f dada, 
como ya se había planteado en la definición de la 
transformación KL. Con este mismo razonamiento 
puedeverse que la matriz de covarianza de laimagen 
de salida puede expresarse como 


Lg) + eE, ASIA 0 


(6.127) 


dado que (1111! - 114". esta última expresión se 
puede reescribir de la siguiente manera 


IS 


(6.128) 


pero la cantidad entre paréntesis cuadrados es 
precisamente la matriz de covarianza K, portanto se 
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que es precisamente la parte medular de la definición 
de la transformada KL. En esta definición se está 
precisamente requiriendo que (1,| seauna matriz 
diagonal para que las bandas de la imagen y estén 
decorrelacionadas. Para que esto sea posible, se 
sabe del cálculo matricial (Arf«en, 1970), que (AJA 
debe ser precisamente la matriz transpuesta de 
los elgonvectores de || siempre y cuando (1) sea 
Una matriz ortogonal, es decir se debe cumplir que 
AA = LObienA* = Al, donde! esla matriz identidad 
y A! es el inverso de A. Este inverso existe 
Generalmente para una imagen mutiespectral 
porque [X;] es cuadrada, simétrica, real y cuyo 
determinante es normalmente diferente de cero. 
Como resultado de estos aspectos (1, — puede 
identificarse como la matriz diagonal de los 
sigenvalores del; portanto 





(6.130 


donde , sonlos eigenvalores mencionados de una 
imagen muhiespectral de y bandas y que esla forma 
explicita de la ecuación (6.115). Hasta ahora la 
transformación KL requiero, desde el punto de vista 
Computacional, de la evaluación de la media de la 
imagen de entrada (ec. 6.122), de la matriz de 
covarianza (ec. 6.124) y de su respectiva 
diagonalización (ec. 6.129). La diagonalización de 
IK;] implica la solución de la ecuación de 
ú'xgenvectores y elgenvalores (ec. 6.115) por medio 
delascuación secular: 
cuando no se requi 
conveniente verificar que exista A”! ya que sinoes así, 
la transformación KL no puede llevarse a cabo, Una 
vez obtenidos los eigenvectores, se construye la 
matriz y de aquí se puede realizarla transformación 
“empleando la relación (6.111), por lo que finalmente 
lo único que se refere explícitamente son los 
sigenvectores asociados a la matriz '*:! , lo demás 
¡son pasos intermedios o de referencia. Además de lo 
anterior, se puede pensar en lo siguiente: como la 
transformación KL es ortogonal, entonces 














(6.191) 


es decir, 





De aquí se puede calcular el porcentaje de la 
varianza »' de las bandas resultantes de la 
transformación KL. Para esto basta dividir cada 
úsigenvalor porel valor de la traza, esto es 





de tal forma que la varianza es un vector de tantos 
elementos como bandas tenga la imagen 
muhiespectral. De acuerdo a lo que se mencionó. 
anteriormente, la correlación entre bandas puedo ser 
muyaka, porloquelos primeros olomontos del vect 

:  tepresentarán el mayor porcentajo de la 
información contenida en la imagen, on 
contrapanida, los últimos elementos 
de representarán el menor porcentaj 
información, es decir prácticamente el ruido de la 
imagen; en otras palabras la transformación KL 
también separa la información del ruido de la imagen 
(Tou y González, 1974) 

Al llevar a cabo la descomposición en 
componentes principales, es muy posible que la 
¡imagen de salida contenga valores negativos que no 
esconveniente manejar nipara ofoctos de despliegue 
ipara transformaciones ulteriores, Debido a esto. 
necesario realizar un reescalami 
dinámico de valores combinado posiblemente con 
una normalización del histograma para darlo a la 
imagen de salida el mejor 


de realce global que tiene esta transformación. Es 
laro, por otro lado, que al aplicar una normalización 
al histograma y un reescalamiento no violamos la 
reversiblidad de la transformación KL, ya que 
aquellas son también reversiblos. Además de esto 
podemos decir que latransformación KL presentalas 
siguientes características: a).- emplea todos los datos 
de entrada, b)- no requiere de parámetros de 
entrada, c)- las tres primeras bandas pueden ser 
eres pera les depa con col 

algoritmo correspondiente es de fácil implantación y 
de rápida ejecución, e)- la normalización 
mencionada puede llevarse a cabo en paralelo a la 
transformación en un mismo algoritmo, f.- en el caso 
de imágenes Landsat, la primera banda de salida es 
proporcional al albedo yla segunda ala primera de la 
transformación Kauth-Thomas, 9).- puede emplearse 
la primera o la segunda banda como representativa 
para realizar alguna otra operación como fitraje o 
realce de bordes y, h).- para efectos de clasificación 
espectral, pueden utilizarse únicamente aquellas 
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bandas que representen la dimensionalidad 
intrínseca de la imagen multiespectral. En el caso 
específico de imágenes Landsat otras aplicaciones 
¡son posibles en las áreas de Geología, Uso del Suelo 
e Hidrología. Puede mencionarse también que la 
descomposición en componentes principales se 
emplea para compresión de imágenes en procesos 
de almacenamiento y envío de éstas. Dado que esta 
transformación reduce la dimensionalidad intrínseca 
delaimagenmutiespectral la representación de ésta 
puede hacerse empleando no todos los 
sigenvectores de la matriz (1, sino Únicamente 
aquellos que comesponden a los mayores 
eigenvalores, es decir aquellos para los cuales el 
porcentaje de varianza sea digamos de al menos 
95%. En muchas ocasiones puede suceder que la 
correlación entre bandas de una imagen 
muliespectral sea tan grande que la primera banda 
de la imagen de salida tenga 95% o más de la 
información total. Debido a esto es entonces posible 
afirmar que esta primera banda contiene los rasgos 
más. relevantes de los patrones presentes en la 
imagen, por este motivo, la descomposición en 
Componentes principales puede verse como una pre- 
clasificación espectral de la imagen original y, por 
sto,también, es que se dice que estatranstormación 
constituye un realce global a la imagen. Por último 
mencionaremos que una aplicación muy interesante 
de esta transtormación , es en el análisis 
mulitemporal de imágenes digitales para efectos de 
evaluación de cambios espaciales ocurridos en el 
tiempo de una O varias clases de patrones. 
Debidamente normalizadas por efectos de 
degradación, ruido y deriva en respuesta de los 











Fig.6.8. Nube de datos de una imagen de dos 
bandas mostrando tres clases separables 
completamente en el eje b,. 








sensores respectivos, o algún otro mecanismo de 
interferencia, las imágenes que contienen el cambio 
mencionado, son registradas unas con otras en un 
sistema de coordenadas común. Una vez hecho esto 
se crea una imagen mutibanda y mulitemporal, 
¡donde las bandas respectivas están compuestas por 
las imágenes adquiridas en diferentes tiempos de la 
escena que se desea estudiar. A partir de esto puede 
plantearse la siguiente hipótesis: habrá una alta 
Correlación entre bandas de la imagen multtemporal 
para aquellas regiones donde no haya habido 
cambio significativo y una correlación baja para las 
z0nas dondelas clases de patrones hayan cambiado 
sustancialmente. Entonces dada la naturaleza de la 
transformación KL, se tiene que el mayor porcentaje 
(ec. 6.133) de la varianza en una imagen 
mulitemporal estará asociado con las clases de 
patrones que presenten poco o ningúncambio, porlo 
que aquellas regiones de cambio localizado 
quedarán diferenciadas o reslzadas en las primeras 
componentes de la transformación KL de la 
mencionada imagen mulitemporal. Cabe mencionar 
también que las bandas de entrada que componen a 
esta imagen pueden, a Su vez, ser las primeras 
componentes de imágenes mutiespectrales 
transformadas por descomposición en componentes 


67.- Análisis Canónico 

La transtormación KL es una transformación de 
similaridad (Arfken, 1970) ya que lleva a cabo un 
cambio de base vectorial que genera un nuevo 
espacio donde setiene una nueva representación del 
operador de covarianza 


y endonde en estanuevarepresentación, el operador 
de covarianza es diagonal con las consecuencias ya 
anotadas. También sabemos que eloperador +) es 
hermitano (real y simétrico), que puede ser 
diagonalizable por el operador unitario ortogonal: ay 
y como elespacio donde estádefinido (14) esfinto 
(de dimensión igual al número de bandas) entonces 
se tiene que sus eigenvalores +, son linealmente 
indepencintas se deironogonaldsynormalzaties 
a la unidad, por lo que generan el espacio 
correspondiente alatranstormación KL. Entonces, en 
esta transformación se está cambiando la 
representación operado IX! auna nueva cada 

1) por medio de la relación arriba anotada, 
Ñaciendó resañar (50. 6:7) que los vectores que 
¡generanlanube de datos de laimagen no cambian en 
su foma intrínseca. En contraparida, en una 
transformación canónica, sicambia larepresentación 
del operador %+! ysí cambia la representación de 
stos vectores (fig. 6.8). El análisis canónico es una 
variante delatransformaciónKL, perocon uncarácter 
más específico pues incluye información acerca del 
comportamiento estadístico de las clases de 
patrones de la imagen muhiespectra y por lo mismo 
se requiere de información a-prioi acerca de la 
misma. En el análisis canónico se tiene una 
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transformación reversible (si se conserva la 
información de las clases) pero no unitaria, en donde 
sereduce estrictamente (no intrínsecamente comoss 
la KL) la dimensionalidad del espacio. donde se 
encuentran representados los datos dela imagen. En 
esta transformación se genera un nuevo espacio 
donde las clases de las cuales 50 proporciona 
información, tienen máxima separabilidad, Se puede 
emplear también el análisis canónico para separar 

¡aclase dada del resto dela imagen, que asuvez se 
ve comola claserestante, Por ejemplo, en lafigura6.8 
+s fácil apreciar el esquema de dos clases de una 
imagen de sólo dos bandas. Las clases a, bycno son 
enteramente separables sobre ninguna de las dos 
bandas originales b, o b,; se requieren ambas para 
producir una eventual separación. Ahora bien, si se 
'e8coge una nueva banda dada por el ojo, (ig. 6.8) 
contal dirección que la máxima separación entre las 
clases a, b y c se proyecte sobre este nuevo ejo, se 
tondrá un método donde se proyecten y separan las 
clases. Cabe mencionar que existen muchas 
ocasiones en que esta separabilidad no se alcanza 
sobre la primera componente de la transformación 
KL, puas la dispersión de las clases y eu posición 
rospectiva on el espacio de muestras puede ser tal 
que aun persista untraslape entre ollas, Dosde luego, 
cuando el número de clases consideradas aumenta 
lasoparabllidad completa lega a sor imposible. Por al 
contrario el análisis canónico funciona muy bien 
¡cuando se desea segmentar una sola clase del resto 
de la imagen, En el marco de lo planteado hasta 
ahora, una transformación canónica: 
sigue 





























donde la connotación vectorial os similar a la que so 
apuntó para la relación (6.111). La matriz. D. roaliza 
la siguiento operación 











e) 
y como consecuencia de (6.135) se tiene 
OS 


6.196) 


En donde Kf esla matriz de covarianza 
promedio de la covarianza de cada clase, es decir 











me 
re 
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(6.137) 


dondaN es el número total do clases consideradas y 

11! 1. esla matriz de covarianza de cada clase 
dada'La matriz — ¡¡!] — representala covarianza 
tre las clases, es'decir es una medida de la 
separación que existe entre allas, Como ya se anotó 
arriba, (1! esla matiz identidad, Hay que remarcar 
que esta fomulación del análisis canónico es 
análogaa la que se hizo para la transformación KL y 
como tal es posible llevar a cabo los pasos que se 
hicioron de la ecuación (6127) a la (5.120) 
introduciendo, desde luego, las condiciones 
específicas planteadas para el anáiis canónico. En 
sl caso de la transformada KL, ora fácil determinar la 
estructura de la matriz |), ya que ésta leva a cabo 
una diagonalización de la matiz de covarianz 
original. No es tan evidente determinar la estructura 
de lamatriz (51, puesto que ésta no diagonal 
la matriz de covarianza de entrada, sino que la 
convierte en una matiz unidad. Para obtenera A 08 
necesario hacer consideraciones explícitas sobr 
Comportamiento estadístico de las clases que se 
desean separar; para esto regresemos ala figura 6.8. 
La idea básica de esta transformación consiste en 
definir un nuevo sistema de coordenadas, sobr 
primar eje del cual se tengan dos elementos: 
Mínima dispersión de las clases consideradas y, la 
máxima separación entro ollas, os decir, se propone 
maximizar las distancias dd, y minimizar x,,x, yx, 
dol ejemplo plantado, Al imponer la condición 
mínima dispersión de ciasos se está requiiendo que 
fticamento óstas sean esféricas, lo que está 
expresado en la relación (6.135). Lina manera de 
roalizar estas dos condiciones es a travós de las 
matices de covarianza de cada clase y la matriz de 
¡covarianza de las clasos como un todo; desde luego 
laprimeraesuna medida dol dispersión delos datos 
para cada grupo o clase considerada y la otra es 
folativa la distancia que oxisto entro los grupos de 
datos. En términos cuantitativos esto se puede 
oscribir como sigue 
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(6.139) 





donde d ss el vector que define el eja del nuevo 
espacio (b, en la figura 6.8) para el cual se tiene la 
ima separación entre clases y la mínima 
dispersión de las mismas. Las cantidades 0; y 
:, son las varianzas, sobre este nuevo ojo, de las 
clases y dela separación entre ellas respectivamente. 
De esta forma, en esta transformación canónica, se 
tiene la máxima separación en el primer eje (b), la 
siguiente separación en importancia en el segundo 
jo (b,) y así sucesivamente. En esto orden de ideas, 
d es el primer sigenvector (aunque después se 
tomará como cualquiera de ellos) y + es portanto 
el primer eigenvalor. La maximización de lo 
axpresado en (6.138) se realiza por métodos 
tradicionales del cálculo diferencial, es decir, se 
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6. CONCLUSIONES 


+ Por primera vez se pudo obtener un mapa forestal de 
escalas 1- 25000 y 1. 50.000, en el cual sa pudieron 
diferenciar ipos y subhipos de bosque, utiizando 
imágenes de radar de alta resolución, (STAR 1 
INTERA), en condiciones del bosque húmedo tropical 
con características de homogeneidad en cuanto a su 
composición florística 

+ Larapidez enla adquisición de los mosaicos de radar, 
STAR 1, (30 días), permitió obtener información 
lorestal oportunamente de acuerdo a la urgente 
nocosidad de evaluar y estudiar la disponibiidad de 
recursos con propóstos de aprovechamiento y 
manejo forestal, lo cual con otro sensor no hubiese 
sido posible en lan corto iempo. 

+ La eficiencia obtenida con el sistema de radar de 
INTERA, para la elaboración de mapas forestales. 
comparada con fotografias abreas de escala 1 
40.000, que soría la escala adecuada para obtener 
resultados similares, es 5.4 veces mayor en tiempo 

+ La posiblidad de obtener mosaicos de radar del 
sistema INTERA en un rango de escaía variable 
(1:50.000 - 1: 200.000). faciitó considerablemente la 
Interpretación. permitiendo correlacionar las 
variables que intervinieron en este proceso 
Igualmente, faciltó un andiess de conjunto y una 
programación detallada delas actidades forestales 
do campo. 

+. La delneación de drenajes pequenos, principalmen- 
te en las áreas planas (lanura aluvial) es clic en el 
sistema de alta resolución (INTERA), debido 
principalmente a la densa cobertura boscosa, 
caractoristica de la zona estudiada 

+ Ladelmilación deuso de latierray de tipos de bosque 
en la región de colmas bajas y altas con radar de 
INTERA, es muy complcada y a veces imposible, 
debido a los contrastes de luminación en las 
pendientes que se encuentran mirando hacia el radar 
(áreas muy claras) y las pendientes opuestas 
(oscuras). 

+ Elrelomodela energía delradar está nfluenciada por 
larugosidad de la superficie vegetal (morfología de la 
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planta), constante dieláctrica, tamaño, densidad del 
folaje y configuración del terreno donde se de: 
sarrolan las comunidades vegetales. Por consi. 
guiente, con ll de identíicar los diferentes tipos de 
vegetación sobre las imágones de radar, el intórprato 
debe tener conocimientos de la reflexión dl la energía 
sobre la vegotación, vasta experiencia en 
totoiterpretación. disponer de fotografías aéreas 
corvencionales o de pequeño formato, las cuales 
serán utilizadas como patrones para clasificar los 
tipos y además conocer la estructura y fisonomia de 
estas. 

Se espera que el mejoramiento de los nuevos 
ssiemas de radar espacial, con mayor resolución y 
“acceso a los usuanos, incremente el uso de estas 
imágenes en los lovantamiontos de los bosques 
húmedos tropicales especialmente para su monitoreo. 
y viglancia 

Las investigaciones adelantadas por ol SIR:B, 
utiizando diterentes ángulos de depresión, han 
aportado conocimientos que permiten comprender 
mejor la relación de la rafiexión de la energía con la 
vegetación y longitud de onda. 

Con las imágenes de radar comparadas en esto 
nico, es posible elaborar mapas de vegetación de 
lipos y subtipos de bosque, los cuales pueden ser 
ullizados para levantamientos exploratorios. (SIR-A), 
de reconocimiento, (WESTINGHOUSE). y semi 
detallados, (STAR-1 INTERA), El uso del radar de 
visión lateral reduce considerablemente el tiempo de 
complacón del mapa y costos en los levantamientos 
de grandes áreas. 

Los sistemas de radar aeroiransportados con 
longitudes de onda menores, permitan la delinaación 
de biomas, ecosistemas, pos y subtipos de bosque 
En el radar de resolución media (Westinghouse). los 
Tipos de bosque no pueden ser diferenciados 
mientras. que algunos sub-pos de vegetación 
pantanosa, en especial los de la llanura aluvial de 
inundación, sí pueden ser determinados. 

Sobre las imágenes del SIR-A, con Imiaciones se 
puede llegar a determinar algunos tipos de 
vegetación menor, poro no tipos de bosque. 





SISTEMAS DE RADAR 
ESPACIALES DE ALTA 
RESOLUCION 





So espera que el Saiáiio. 
ERS-1, de próximo lanza 
miento por la Agencia 
Espacial Europea (ESA), 
aporte importantos posibi 
lidades alestudo de zonas 
de vegetación tropicales, 
dadas su alía rosolución 
“espacial (30x30 ms) y su 
capacidad de loma de 
'imágonos en condiciones 
de alta nubosidad. 
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esla poblacióntotaldelas clases consideradas y Mes 
elnúmero de éstas. Porotrolado, auncuandosetiene 
a la ecuación (5.125), como referencia, conviene 
escribir explíctamentola expresión correspondiente 
para 1551 


hi A 
UL lr 004 06, 


(scordando que se obienade a expresión 6122) 
'esto queda completo el esquema para ejecutarla 
transformación canónica aquí planteada, hacemos 
ónfasio, desde luego, en que es necesario conocerlas. 
clases involucradas, es decir, hay que contar con los 
valores y número de los pixels que las componen. 
Esto se obtiene claramente a parti de la imagen 
original y de regiones rodalizacas de ésta donde se 
pionse que existe la clase o clases de interés; este 
trabajo puede hacerse también a partir del anáíisis de 
cúmulos (clustering) de la imagen mutiespectal. Es 
necesario mencionar también que los oigenvalores d 
de esta transformación no son en general 
rtonormales pero pueden hacerse así empleando el 
método de Gram-Schmich (Arikon, 1970). 


6.8.- Transformación Kauth-Thomas 


Un ejemplo relevante de una transformación 
canónica aplicada específicamente 
'multiespectrales de la serie Landsat, con aplicación 
particular a evaluación de cobertura vegetal agrícola, 
+s la que se conoce como Transformación Kauttr 
Thomas (KT), diseñada por estos investigadores 
(Kauth y Thomas, 1976) para los primeros Landsat y 
extendida posteriormente al TM de los últimos 
satóltes de este serio (Crist. 1984). 

Con el objeto de ilustrar el diseño de la 
transformación KT, roficámonos a la figura 69, en 
donde se muestra esquemáticamente la nube de 
datos que se relaciona con un cultivo agrícola con 
plantas en diferentes etapas de crecimiento. Para 
poder mostrar correctamente la evolución en el 
tiempo de un tipo de cultivo a partir de una sola 
imagen, es, desde luego, necesario suponer válido el 
teorema ergódico para la imagen muliespectral 
Correspondiente. En la figura6.Sa, sedalagráfica, en 
elespacio generado porlas tres primeras bandas (b,, 
b,, b,) de una imagen Landsat (MSS), de una 
¡subregión de ésta donde se sabe que existe un cultivo 
determinado. En el desarrollo de un cultivo podemos 
pensar en el siguiente proceso: en primer lugar a 


través de extensivas medidas sobre muestras de 
suelos (Cond. 1970), ha quedado demostrado que 
los respectivos datos de refectancia, quedan 
aproximadamente sobre una línea recta cuando son 
graficados en un espacio formado por las bandas 2y 
3 del MSS de la serio Landsat (Kg. 6.90), Sobre esta 
línea, los suelos oscuros quedan hacia el origen de 
coordenadas y los suelos brilantes sobre el otro 
sxtremo de la línea; los diferentes tipos de suelos 
Quedarán desde luego sobre puntos intermectios, Al 
observar la evolución espectral de un determinado 
Cultivo, se tendrá en primera instancia la mencionada 
línea de suelos al momento en el que las plantas 
respectivas apenas comienzan a aparecer sobre la 
superficie de la tierra. A medida que las plantas 
crecen, la respuesta espectral de éstas sa combina 
proporcionalmente con la del suelo donde sa 
encuentran plantadas, lasombra que éstas proyectan 
sobre el suelo y su propia respuesta irán dominando 
sobre la de éste, de tal manera que lo que antes era 
suelo brilante aparecerá ahora como una mezcla 
másoscura(fgs. 6 Say8.9b). Paralos suelos oscuros 
la situación se mantendrá aproximadamente estable 
por un tiempo en cuanto a la mezcla de respuestas: 
uelo-plantas-sombras. Cuando las plantas emergen 
decididamente sobre la superficie terrestre. SU 
respuesta espectral empieza a dominar sobre los 
tros factores (suelo y sombras) puesto que su ároa 
Yoliar cubrirá paulatinamente al CIV correspondient 
de esta manera las respuestas combinadas de 
plantas sobre suelos oscuros y brilantes se 
aproximarán a una sola en el punto de madurez dol 
Cultivo (figs. 6.3a y 6.90). De aquí que la evolución 
arriba mencionada, describa una gráfica que se 
asemeja a la de un penacho. Después de alcanzar la 
madurez, las plantas entran en la etapa do 
senescencia y empiezan a tomarse amarilas, 
describiendo, su respuesta espectral, na trayectoria 
Curva hacia la línea de suelos; esta trayectoria tiene 
una ciena dispersión (fg. 6.9a) dependiendo de 
ciertas diferencias que pudieran tenor las plantas del 
¿uttivo, o bien si existen varios cultivos semejantes en 
la misma área de plantación. La variación on 
respuesta, producto de la senescencia, agrega a la 
forma del penacho una especie de cola o bora y de 
aquí el nombre que a veces se le impone a la 
transformación KT: gorro de bora, Síbien el proceso 
de emergencia y madurez del cultivo describe la 
forma de una especie de triángulo aplanado, 
básicamente sobre el plano de las bandas b, -b,, la 
tapa de senescencia se sale de este plano (f.6'9c) 
por lo que se requiere finalmente de un espacio 
tridimensional formado por las bandas b,-b,-b,; esto 
Quiere decir que la transtormación KT genera un 
"nuevo espacio de tres bandas que es suficiente para 
realzar convenientemente el proceso evolutivo de un 
Culivo determinado. Cabe aclarar que lo que se ha 
¡denominado plano de los suelos (figs. 62ay6.90) en 
realidad es una nube de datos formando un elipsoide 
muy aplanado y bastante elongado. Dado que en 
General puede afirmarse que cualquier especio 
vegetal tiene un ciclo de emergencia, madurez y 
senescencia similar al ariba descrito, es posible 
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Fig. 6.9. Descripción esquemática del proceso de crecimiento de un cultivo a pantir de datos Landsat. 











afirmar que la transtormación KT es en realidad un 
realce para cualquier cobertura vegetal Es 
importante. puntualizar sin embargo, que esta 
transformación requiere de datos de entrada 
obtenidos de aquella subregión que contenga la 
comunidad vegetal de interés, por lo que el realce 
seráenfático conrespecto ala comunidad estudiada. 
Para finalizar este párrafo diremos que los detalles 
físicos sobre el proceso evolutivo mencionado 
pueden encontrarse directamente en el artículo 
¿original de Kauth y Thomas (1976). 
LatranstormaciónKT, al igualque la KL y el anáíisis 
canónico, pueden verse como una pre-clasificación 
de una O varias clases de objetos o rasgos de la 
escena, de aquí que estas transformaciones se 
emplean a veces como datos de entrada para un 
dlasificador espectral convencional, con reducción 
sen costo de computación ya que la dimensionalidad 
de la imagen original siempre se reduce. Sin 
embargo, para clasificaciones absolutas o para 
comparaciones ya sea mutitemporalo mutiespacial, 
es necesario incluir varios aspectos extemos a la 


evaluación de la cobertura vegetal. Estos aspectos 
extemos son básicamente de dos tipos: condiciones 
atmosféricas y geometría de iluminación. En cuantoal 
primer efecto, es claro que componentes como 
neblina, bruma, aerosoles, nubes, vapor de aguay la 
distribución del albedo en la zona de interés, afectan 
sensiblemente la medida de ragiancia hecha por el 
sensor sateltario y en consecuencia la respuesta 
espectral observada no corresponde directamente a 
los objetos sensados y además puede variar de lugar 
alugar y de untiempo a otro. Como ya se mencionó 
con anterioridad, las sombras que proyectan las 
plantas en su etapa de crecimiento y maduración, 
ejercen una influencia en la variación de la respuesta 
espectral observada. Además de esto, a medida que. 
cambian los ángulos de vista e iluminación, la 
cantidad de sombra que proyecta el cultivo cambia 
también; de aquí que la refectancia total pueda ser 
más clara u oscura de lo que en realidad es, Es por 
esto que la geometría de iluminación es un factor 
importante en la observación de cultivos por medio 
dela transformada KT. Debido a estos dos factores 
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externos, es importante llevar a cabo 
normalizaciones por efectos atmosféricos (Lira y 
Oliver, 1983) y por cambios en iluminación antes de 
aplicar la transformación KT, esto como ya se 
mencionó, en el caso de que se deseen medidas 
absolutas de cobertura vegetal. 

Veamos ahora los aspectos computacionales de 
esta transformación; ésta se hace de la siguiente 
manera. 





dondela notación vectorialtiene el sentido explicado 
anteriomente, fr) es la imagen multiespectral de 
entrada, g(u) la de salida, R! es la transpuesta de la 
matriz de transformación y r, es un vector arbitrario 
para evitar valores negativos en la imagen de salida, 
aunque si bien este vector puede evitar valores 
menores a cero, también puede ocasionar algunos 
Puntos de saturación en la imagen de salida, por lo 
'Quees másrecomendable realizar un reescalamiento 
lineal del histograma de g(u). La matriz R es unitaria y 
se construye de la siguiente manera. 

El primer vector columna R de ¡| se escoge de 
tal manera que su dirección esté alo largo del eje 
'mayordelelipsoide queformala nube de datos de los 
suelos dela subregión o subregiones seleccionadas, 
recordando que estos datos se obtienen a parir de la 
imagenmisma. Enla práctica R, se obtiene a parir de 
la descomposición en componentes. principales 
aplicada al elipscide de datos y escogiendo su 
dirección en coincidencia con la de la componente 
principal. 


Elsegundo vectorcolumna R, de (21 ,se escoge 
perpendicular a R, y apuntado en la dirección 
principal de maduración del grupo vegetal 
Seleccionado. Esto se hace a partir de proyecciones 
de la nube de datos como la que se presenta en la 
figura 6.9b, en donde se aprecia dicha dirección; la 
descomposición en componentes principales 
también juega un papel aquí. 








Eltercer vector columna R, de _ (1), se escoge 
desde luego perpendicular a, y R, y apuntando en 
la dirección de senescencia (fg. 69c) de la 


Comunidad vegetal de interés. Nótese que en las 
subregiones seleccionadas de la imagen para servir 
¡Como datos de entrada, debe haber plantas en todas 
tapas de su evolución y una variedad suficiente de 
suelos, de otra manera los vectores R, R, y R, no 
pueden ser estimados con precisión. 


El cuarto vector columna R, de 1! , se escoge 
obviamente perpendicular a os tres anteriores. 


Una vez evaluados estos cuatro vectores, para 
imágenes de cuatro bandas del MSS de Landsat, se 
aplica el método de ortonormaización de Gram- 
Schmidt con el objeto de conformar a (| como una 
matriz unitaria yortogonal. De acuerdo alo quese dijo 











Fig.6.10. Barrido de una escenaa incrementos 
angulares iguales: el CIV del centro es más 
pequeño que elde los extremos. 





anteriormente, esta transformación fue extendida a 
imágenes de siete bandas del TM de Landsat (Crist, 
1984; Crist y Cicone, 1984), pero la idea es 
básicamente la misma que la aquí expuesta. Los 
vecoresh,. Re Ry, generanun nuevo espacio de 
dimensionalidad igual a tres, ya que en la cuarta 
banda de salida prácticamente no hay información 
acercadelacobertura vegetal en general, Ala primera 
banda de la transformación KT, dada por R,, se le 
lama brilantez y es proporcional al albedo. A la 
segunda banda, dada por R,, se le llama verdura 
(greeness) y da una medida del vigor de la 
vegetación. A la tercera banda, dada por R,, se le 
llama amarilez (yellowness) y da una medida de la 
senescencia del cultivo. A la cuarta banda, dada por 
A, se le llama ausencia (non-such) y tiene este 
nombre porque carece de información relevante 
acerca del culivo, Entonces el espacio generado por 
los primerostres vectores puede ser empleado como 
un medio para evaluar el estado de uno o varios 
Cultivos o para estudiar la cobertura vegetal de una 
región dada. 


8.9.- Transformaciones Geométricas 


Enel capítulo anterior se discutió la digitización de 
una escena empleando un formalismo. vectorial 
conocido como lattica de muestreo, La imagen 
resultante queda por tanto implícitamente referida al 
sistema de coordenadas definido por la latica 
"empleada. Existen, sin embargo, factores externos al 
proceso de muestreo que pueden hacer que el 
sistema de coordenadas al cual realmente se refiere 
la imagen no sea precisamente el de la latice, 
Algunos de estos factores son los siguientes: la. 
"escena no es un plano sino más bien una superficie 
cualquiera, la escena es un plano pero el sensor la 
barre en incrementos angulares constantes (fig. 
8.10), poro que los CIV extremos son más grandes 
quelos del centro, existe un movimiento relativo de la 
escena en la dirección de barrido del sensor, la 
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plataforma del sensor no es perfectamente paralela a 
la escona, existen aberraciones de consideración en 
el subsistema Óptico del sensor, el CIV no es 
cuadrado sino que tiene alguna otra forma: 
rectangular, cónica o romboide, existe algún efecto 
de dispersión no-ineal que distorsiona la imagen y. 
'en general, cualquier interferencia en el sensor y la 
escena o alguna no-inealidad en el funcionamiento 
de aquél. Además de la posible corrección de estos 
efectos, es factible que se requiera proyectar la 
imagen a algún otro sistema de coordenadas donde 
resulte más conveniente su análisis, Debido a lo 
anterior, es necesario establecer un procedimiento 
paratransformar el sistema de coordenadas alcual se 
roflore originalmente la imagen; a este tipo de 
transformaciones se les conoce como geomélricas 
porque precisamente cambian la geometría de la 
Imagn. Algunas de estas transformaciones pueden 
roalizarse modificando la geometría de captura de la 
escena. introduciendo un hardware" especializado 

nel subsistema de formación yregistro de laimagen. 
La mayoría, sin embargo, deben hacerse a-posteriori 
por medio de un modelo matámatico bien definido, o 
empleando polinomios de un cierto grado en 
combinación con lo que se conoce como puntos de 
control. Es necesario aclarar que en la literatura se 
habla en ocasiones de correcciones geométricas 
para compensar por uno o varios de los factores 
mencionados. Para hablar de correcciones es 
necesario suponer que el sistema de coordenadas al 
cual so quiere mapear la imagen, es precisamente el 
correcto. Esto sólo puede hacerse si se introducen 
elementos de valor en estas transformaciones, es 
decir setione que asumir un sistema de coordenadas 
correcto porque de hecho es en éste donde 
asignamos un mayor valor o utilidad a la información 
contenidasn la imagen, o bien, 
'elnuevo sistema de coordenadas la imagen tiene un 
aspecto mejor para propósitos de inspección visual. 
Con esto en realidad se admite que la imagen como 
un todo o ciertas clases de 
realizadas on el nuevo sistema de coordenadas, En 
lortas ocasiones se tienen sistomas de coordenadas 
universalmente aceptados por lo que es conveniente 
roforirla imagen a uno de óstos ya que a panirde aquí 
so pueden hacer mediciones o comparaciones 
cuantitativas de varlos aspectos geométricos de la 

Por último, diremos en esta discusión que 

cualquier transformación geométrica debe hacerse 
banda por banda de la imagen muhiespectral. 

Una transformación geométrica se expresa 
matemáticamente de la siguiente manera 


























(6.152) 


donder = (|) es elvectorde posición de los pixels de 
laimagenen elsistema de coordenadas original, w = 
(uy) es el vector de posición de los pixels de la 
imagen resutante y 17) es la matriz o modelo 
matemático de la transformación geométrica 





(Richards, 1986). A veces es conveniente escribir la 
versión escalar de la ecuación (6.152) 


EI 
ao 
(6.150) 


Hay ocasiones en que el modelo exacto de la 
formación geométrica no se conoce y entonces 
idecuado emplear aversión escalar aproximando 
los operadores /, y /, por un polinomio que 
góneralmente se trunca para incluir sólo los términos 
de orden dos. De esta manera se genera una 
aproximación para u y v en función de | y k. 
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(6.154) 


donde las a, y by son constantes por determinar en 
cada caso de transformación geométrica. Este 
constantes se determinan a partir 
control escogidos éstos de tal 
minimice el error entre el valor verd. 
(us) y ol estimado ¿%.;) .Los puntos de control son 
Puntos observados o identificados sobre la imagen 
original y que se sabe por algún método, que 
corresponden a puntos en el espacio de salida y 
cuyas coordenadas en este nuevo sistema son bien 
conocidas, 

Los puntos de control se escogen lo más 
uniformemente distribuidos sobre la imagen para 












ostimación polinomial. Si se tienen M puntos de 
ontrollacorrespondencia.n el espacio dolaimage 
(fo. 6.11) se expresa en forma vectorial como sigue 





A] 
(6.155) 


Nótese que hay una correspondencia directa entre. 
los puntos do control (|, k,) y ostos que son los 
observados, es decir los puntos de control se 
idontican en la imagen original y en la rejila: de 
referencia (figs. 6.11 y 6:12). De manera similar se 
pueden escribir las constantes involucradas en la 
aproximación polinomial 








EA] 
ba [besbio o ++ ba) 
(6.156) 


Al emplear esta notación, es fácil ver que el error 
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Fig. 6.11. Transformación geométrica de una imagen (izquierda) por medio de puntos de control w 
observados (cruces) sobre una rejilla de referencia (derecha). 








'medio al cuadrado tiene la siguiente forma 


End AM) A) 0 (AD G  a 


(6.157) 
suponiendo claro está que el error sobre una 
dirección del espacio de coordenadas es 


independiente de la otra. La matriz A está definida a 
continuación 


Dai bd 
Ml Lia 


DA 


(6.158) 


Esdecir,lamatriz. (11 se construye a partir de los. 
Pantos de cotrl (Je) Kdentcados en arejla de 
referencia (fgs.6.11y8:12). Además de esto, es bien 
¡conocido de la teoría de operadores lineales (Pran, 
1978) que el error mencionado es mínimo si 


»- ala 


»oale 


donde 11!) es elinverso generalizado (Rao, 1971) 
de (4). Porotro lado, también es fácil er (Pratt, 1978) 
que si el número de puntos de control es mayor al 
número de coeficientes del polinomio empleado, 
entonces 


ah catala 


(6:19) 


Los puntos de control (ec. 6.155) se introducen en 
la aproximación polinomial (ec. 6.154) para obtener 
los coeficientes (a, -*,», si el número de puntos de 
controles mayoralñúmero de coeficientes setieneun 
sistema de ecuaciones sobredeterminado y los 
coeficientes se pueden calcular Óptimamente por 
algún método de mínimos cuadrados. Para un 
conjunto dado de puntos de controlse puede evaluar 
elerror correspondiente (ec. 6.157) y compararlo con 
sel tamaño asociado a la mínima celda de la rela de 
referencia (Kg. 6.11); si el error es razonable para un 
Problemadado entonces se continúa con el proceso. 
Enla aplicación de una transformación geométrica 
¡generalmente se tiene ala mano un número finito de 
puntos de control de tal manera que no es posible 
encontrar siempre un subconjunto de éstos que 
minimice el error estimado. Lo que sí se alcanza casi 
siempre es el determinar un juego de puntos de 
control más o menos uniformemente distribuidos 
(Schilen, 1979; Ort 1981) sobrelarejila dereferencia 
'queaseguren unnivelrazonable de error, De aquí que 
noserecomiende nunca elescoger puntos de control 
que estén preferentemente localizados enuna región 
O, en general, que no tengan una distribución 
homogénea e isotrópica sobre la rela mencionada, 
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Al determinar los coeficientes (a,b) la 
aproximación polinomial queda entonces lista para 
emplearse en latransformación geométrica deseada. 
Para esto es necesario realizar un remuestreo del 
espacio de la imagen original al de la imagen 
corregida; el procedimiento es como sigue: la rejila 
dereferencia debe tener un tamaño de celdilaigual al 
tamaño asociado a los pixels de la imagen original, 
esta rejila está referida desde luego a un sistema de 
coordenadas prefado, si el tamaño de los pixels 
originales varía, la celda será asociada con el mayor 
de éstos. Comenzando entonces con la celdilla del 
espacio de salida, cuyas coordenadas en la latice 
original son u = 1 yv = 1, es decir, aquel valor de y 
k para los cuales se tengan precisamente las 
coordenadas del origen de la imagen de la cual se 
parto. Tomando como referencia a este punto inicial 
(k) y deacuerdo alarejilade referencia predefinida, 
86 incrementan secuencialmente | y k, se evalúa la ú 
y Ocorrespondientes, las que no necesariamente se 
felacionarán a una posición exacta (fig. 6.12) de un 
punto de la latico original, porlo que será necesario 
realizar una interpolación a partir de los valores delos 
pixels vecinos. Cuando para un juego de valores y k 
se tenga que las ¡ y 5 correspondientes exceden 
“completamente los límites de la imagen original, 
entonces se asignaun valor cero. pixelde salidayse 
ajustan los valores de j y k para continuar con la 


siguiente línea o columna de la rejila referencia, 
según se haya determinado en el proceso de 
romuestreo. Es claro que de esta manera los límites 
dela imagen de salida serán en general no paralelos 
alos ejos del sistema de coordenadas de larejila de 
roforencia, por este motivo es necesario rellenar con 
ceros hasta tener límites paralelos, Esto además, es 
necesario hacerlo, puesto que computacionalmente 
es más fácil manejar archivos de forma asociada 
rectangular y así visualmente, la imagen de salida 
rosata mejor en una pantalla de despliegue. 

La construcción, en la rejila ideal, de la imagen 
transformada, implica evaluar las coordenadas i y 7 
del pixel que corresponde a un par ()) dado, sin 
'smbargo, las posiciones estimadas (1.%), a partir do 
la aproximación polinomial, no necesariamente 
coinciden (fig. 6.12) con las posiciones de los pixels 
de la imagen original. Debido a esto es necesario 
realizar una interpolación para calcular el valor del 
pixel de salida: aquel con coordenadas (5,5). La 
¡evaluación de este valor puede ser tan sencilla como 
asignarle el que tiene el pixel más cercano en la 
imagen original. Se puede también llevar a cabo una 
interpolaciónlineal a partir de los pixels más cercanos 
(vecinos 4 u 8) al pixel de salida, o bien lo que se 
conoce como convolución cúbica. Podemos decir 
ahora que una transformación geométrica se 
compone de dos operaciones básicas: a) una 
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Fig.6.12, Mapeo de los pixels cuyas posiciones son estimadas porla transformación geométrica, el valor 
del pixel de salida se ostima por interpolación de los valores de los pixel vecinos. 
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Iattice de la imagen original 









Fig.6.13.- interpolación bilineal para obtener el valor del pixel cuya posición estimada es (4,v). 


páxol de solida 





onición estimado (5,5) 








transformación espacial que define elrearreglo de los 
Pixels de laimagon original yb) una interpolación que 
Ínvolucra los valores de los pixels originales con el 
Objeto de generarlos valores delos pixels de salidaen 
únarejila discreta dela imagen resultante la parte (a) 
ha sido ya discutida en esta sección. En cuanto a la 
interpolación se puede escribir, en general, como 
sigue 


donde R es la función de interpolación y ';; es una 
subregión finita de la imagen original, es decir, para 
evaluar el valor de un pixel de salida se toman 
solamente Jo les ca una veciciad feta que 
contenga al pixel de coordenadas (;.; 

proceeds de" ineciación que” meracan' una 
discusión en detalle son/os lineales y de convolución 
úbica; el del más próximo vecino estan sencilo que 
norequiere de explicación adicional, además de que 
en ciertas ocasiones introduce artefactos por lo que 
no es muy recomendable su empleo. En cuanto a la 
interpolación lineal, que en procesamiento de 
imágenes se conoce más bien comointerpolación bi- 
lineal, fijémonos (fig. 6.13) en una vecindad de forma 
¡cuadrada alrededor del pixel cuya posición (3,;)) ha 
sido estimada. Entonces, a partir de los pixels más 


Cercanos, se obtienen los interpolantos 





dia? vol) 4 Lluel 081) 9 (La t10(U,v81) 


(6.161, 


donde se ha supuesto que de pixel a pixel en la 





parti 
u, vyson respectivamente u” y v. El paso final ahora 
esinterpolarlinealmente entre(+',y a(u',v+1) para 
obtener 








RULO = Aaa v 0) (ea. e 64! va) 


(6.162) 


e deci, QUA es el valor interpolado del pixel de 
a pÍAÍióN es 0/0 on tica tal Enla 
rain cs gls ra vera yo 
'enlainterpolación!ineal donde sólo se utiizanlos 
duao plele más próximos (vecindad tio 4) 
Considérense ahora los 18 vecinos más cercanos a 
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entonces sobre las cuatro líneas (fig. 6:14) de 
cuatro pixels que se encuentran en esta vecindad, se 
obtienen cuatro interpolantes empleando polinomios 
cúbicos (de orden tres) cuyo resultado son cuatro 
valores intermedios para las posiciones marcadas 
conxenla figura 6.14. Sobre estas cuatro posiciones 
intermedias se aplica de nuevo una interpolación 
cúbica para obtener finalmente el valor del pixel 
correspondiente a la posición — (5,:). El álgebra 
involucrada en estas, sones es bastante 
elaborada (Shlien, 1979; Simon, 1975) porlo quesólo 
se presentan los resultados aquí. Elinterpolante para 
la primera Ínea, suponiendo que el pixel de la 
izquierda superior (fg.6-14) es precisamente el K(u.w), 
está dado por 








Los otros tres interpolantes. 
y 201043) se obtienen directamente a panir de 
(6.163) variando simplemente elíndice v para las tres 
líneas respectivas, es decir hay que substituir en lugar 
le y, los valores v+1, v+2 y v+3. Con estos 
interpolantes se obtiene finalmente el valor del pixel 


Cu" .vo1), 00 








de salida correspondiente a la posición estimada 
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Tantola interpolación bilineal como la convolución 
cúbicason, a finde cuentas, promedios ponderados 
sobre una ventana finita, esto introduce un fitro de 
pasabajas (Lira, 1989) que puede degradar laimagen 
resultante haciéndola aparecer como en un efecto de 
fuera de foco. Es conveniente entonces aplicar una 
restauración (Lira, Oliver, 1983) para recobrar la 
apariencia de nitidez de la imagen original. 

La aproximación polinomial introducida en esta 
sección incluye polinomios de orden dos 
únicamento, ya que la experiencia ha demostrado 
que éste es un buen compromiso entre la precisión 
requerida y el costo computacional. Un polinomio de 
grado uno es en realidad una transformación lineal, 
insuficiente para introducir las modificaciones 
espaciales que generalmente requiere una 
transformación geométrica. Un polinomio de orden 
tres o más puede ser útil para tomar en cuenta 
modificaciones de orden menor como son las 





lattáce de la imagen original 










Fig. 6.14. Interpolación por convolución cúbica, primero sobre las cuatro filas de pixels vecinos al 
estimado y después sobre la línea marcada con cruces. 


pixel de salida 


FÍN interpolado 


ponición estimada (4.5) 
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poqueñas desviaciones a la geometría ideal del 
sensor con respecto a la escena, Un buen cuidado 
debo tenerse sin embargo, pues sino se introducen 
suficientes puntos de control, un polinomio de orden 
reso mayor, generalmente introduce inestabilidades 
en el proceso de transformación geométrica que 
necesariamente —acarrearán deformaciones 
indeseables en la imagen de salida. Con el objeto de 
ver con mayor detalle estos aspectos acerca de 
transformación polinomial, escribamos nuevamente 
las ecuaciones (6.154) 





Do da Ide bale e gba bal 





108 dela derecha, es decir, (4, b,)ropresenta 
un desplazamiento de la imagen transformada con 
respecto ala original, los términos segundo y tercero 
sorefioran a cambios de escala, los tórminos tercero, 
cuarto y quinto son rotación y redistribución espacial 
do laimagen transformada en relación a laoriginal, Si 








introducen en el polinomio, factores 
entonces se estarán tomando en cuenta 
rodistribuciones espaciales con más detalle para 





efectos menores de la geometría involucrada on la 
producción de la imagen. En resumidas cuentas, 
desplazamiento, cambio de escala, rotación y 
redistribución espacial delos pb imagen, son 
los factores necesarios en cualquier transformación 
geométrica, En el caso de imágenes de satólto 
(Landsato deradar), cuando los puntos de control se 
obtienen de un mapa para una proyección 
cartográfica dada, la transformación goomélrica 
toma el nombre particular de geocodificació 
manera que cada pixel de la imagon se refiere en 
tórminos de coordenadas geográficas. Sisetiene una 
distribución bi-dimensional de datos obtenidos por 
unsensor cualquiera y se relacionan a un mapa por al 
método discutido aquí, entonces el proceso se 
conoce como georeferenciación; así un 
levantamiento gravimétrico o magnetométrico puede 
aciquiri validez cartográfica, lo que es muy úti para 
Integrar datos geofísicos e imágenes de satólite a un 
sistema de información geográfica. El término de 
registración es utlizado cuando dos imágenes de la 
misma escena, pero adquiridas en condiciones 
diferentes, se relacionan al mismo sistema de 
coordenadas. En la registración, las imágenes 
sujetas a transformación, pueden geocodiíficarse 
todas a un mismo mapa con los mismos puntos de 
control, o bien se derivan puntos de control apantir de 
unade ellas y setransforman el resto de lasimágenes 
¡conbase en estos puntos. La registración permite por 
tanto la comparación directa de imágenes de la 
misma escena que hayan sido generadas bajo 
diferentes condiciones de vista de ésta, 














6:10.- Transformaciones Quasi-lineales 


Se han discutido hasta ahora una serio de 
transformaciones reversibles con diversos grados de 
simetría, teniendo las de mayor grado como las 
unitarias hasta las de menor grado como las 
modificaciones. trigonométricas “al. histograma. 
Algunas de las transformaciones requi 
parámetros de entrada, otras no; pero todas con el 
objetivocomún de mejorarla apariencia general dela 
Imagen y sin discriminar a-prior un cierto contenido 
de patrones, esto es así aun en el caso de los 
operadores diferenciales en donde los patrones 
espaciales de alta frecuencia son mayorment 
enfatizados que los demás. Desde luego 

mejoramiento de una imagen está siempre 

















"enmarcado en elcontexto de un problema y tiene que 
ver mucho con juicios subjetivos de valor y calidad 
que.l investigador asigne a su imagentransformada, 
La gran mayoría de las transformaciones estudiadas 





operaciones trigonométricas las únicas que no lo 
son, aunque por su naturaleza son tambión 
irreversibles, No existo un formalismo goneral para 
transformaciones no-iineales como lo hay para las 
que sí lo son; no obstante en ciertas ocasiones 
algunos problemas no-inealos pueden. quasi. 
Iinoalizarse, En el análisis do imágones podemos 
encontrar una clase de problemas no-inoales que 
ropresentan situaciones exporimentalos de una 
importancia que valo la pena su atención, Dos 
enfoques a esto son posibles; uno: 












sde lo que 
originalmente es una transformación no-ineal. 

Los términos linealidad y superposición se aplican 
comúnmente sólo a sistemas linoalos aditivos; sin 
embargo estos conceptos pueden tener una 
connotación más amplia. Oppenheim y 
colaboradores (Oppenheim, 1967; Oppenheim at. 
al, 1968) introdujoron el concepto de superposición 
generalizada, como una manera de extender la 
noción deinealidadacitiva aclasos más goneralos de 
sistemas de los que hasta ahora so han estudiado. 
Para introducir esta generalización, tómenso dos 
imágenes f, y f,, combinadas de acuerdo a un 
operador / que por el momento se pensará como 
indefñnido, esta combinación se supone que produce 
una imagen de salida y 





BUY) Eb) 0 1661) 
(6.165) 


Entonces /1:1  esuna operación Iineal general 
si posee las siguientes propiedades 


Est) 0 Eo(40) = 010.7) 0.0 165,1) 


(6.166) 
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Es decir, que la línealiza la 
combinación entre las dos imágenes , y 1, Desde 
luego esta operación debe de estar detal 


manera que produzca el resultado de las dos últimas 
relaciones. Oppenheim también demostró que si el 
operador O. representa adición en un espacio 
vectorial y si los dos puntos denotan multiplicación 
escalar, entonces el operador “1'] puede ser 
descompuesto en una cascada de operaciones que 
colectivamente son llamadas fitro homomórfico. En 
ste caso las ecuaciones (5.166) y (6.167) asumen la 
siguiente forma 


AN O ANT AN 


dlls LG) MOL) 
(6.169) 


Dela ecuación (6.168) se aprecia que enla primera 
etapa de dicha cascada el operador O es convertido 
sn adición algebraica. En la segunda etapa se tiene 
quo el sistema original, O f, ha sido transformado a 
un sistema lineal acktivo ordinario. En la tercera y 
Úlima etapa de la cascada se aplica el inverso de 
para obtener 


GN 0 
(6.17 


desde luego se supone que existe elinverso de / y 
que —-1% - ¡ , donde: esla operación identidad. 
Ahora bien un sistema generalizado se representa 
esquemáticamente como sigue 








tn] + 0 





donde f, representa la imagen o las imágenes de 
entrada dl sistema, que después de sertransformada 
producenla imagen de salia ,.Elfitro homomórfico 
de dicho sistema generalizado también se visualiza. 
sn forma de diagrama de bloques 


Linear 
aéltico 





ss > [01] 








Elfitrohomomórfico convierte entoncesal sistema 
¿generalizado un uno lineal aditivo, en esta etapa se 
realizan las Operaciones sobre las imágenes 
involucradas de acuerdo al formalismo desarrollado 
en el presente trabajo y posteriormente se regresa 
larepresentación generalizada por medio del inverso 

47, parafinalmente obtener las imágenes de salida, 

Elfitrohomomórfico es, portanto, una herramienta 
para enviar un sistema generalizado a un espacio 
¡donde se puede aplicar la teoría general establecida 
para sistemas lineales aditivos. Es por esto que un 
fitro homomórfico linealiza todos aquellos sistemas 
que en realidad son quaskineales, Es claro que no 
siempre es posible esta linealización ya que no 
úslempres factible encontrara / para que realico esta 
tarea, o bien el inverso respectivo no existe, a menos 
que se trabaje con un pseudoinverso de // (Pratt, 
1978). 

'Uno delos ejemplos más típicos y más socorridos 
del ftraje homomórfico es aquel que se aplica a lo: 
sistemas multiplicativos. Supongamos entonces que 
la imagen olx) se penerada por a mulipicacón 

escalar (pixel a pixel) de las imágenes de entrada 
1,0<y) y Ey), por lo que se tiene 








8007) 0000 0 ln yb) 


(61m) 


En este caso particular es fácil ver que sl tomamos 
/K1 « torl:]_ resulta en la aplicación del logaritmo 
a ambos lados de esta ecuación, por lo qu 
produce un sistema lineal aditivo. 





Lona) 1 dol £a009)) + done] 


(670 


Es decir, en esta etapa, los logaritmos de las 
imágenes de entradason combinados aditivamente y 
cualquier operación al sistema deberá hacerse sobre 

esta baso, Después de haber ejecutado cualquier 





a espacio de jo 
oil SEX, EN un ciagrama de bloques, estos 
aspectos se representan como sigue 





Lineal 
altivo 


Ln + ss 





1o8(-1 




















IAS 


dondeen estecaso, , esla forma esquemática de las 
imágenes de entradaylaimagen de salida!, esla que 
resuñe de todo el proceso que se aplique al sistema 
lineal aditivo. En cuanto a la operación del producto 
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escalar para sistemas muliplicativos, se tiene 
800) ke 10 


donde tes una imagen de entrada cualquiera y k una 
constante real. Al aplicar el logaritmo a ambos lados 
de (6.173) resulta en 


Lonla(2.9)] = klogl£Gx.1)1 


Esta aplicación del fitrajo homomórfico a sistemas 
mukiplicativos es de paricular importancia en las 
técnicas de restauración espacial de imágenes. Este 
tipode restauración hatenido gran éxito paracorregir 
imágones que han sido adquiridas en condiciones de 
iluminación dispareja (La Brecque, 1988), los detalles 
de esto se verán en el capitulo sobre restauración de 
imágenes y se discutirán con ampitud dada la 
importancia que tienen en la comprensión de los 
mecanismos de la psicovisión humana. 

Dosclo el punto de vista físico, un sistema lineal se 
caracteriza con base en un operador que representa 
larespuesta del sistema ala función de impulso, detal 
manera que la imagen de salida os la superposición 
de los impulsos de entrada que pasan a travós del 
sistema, el cual los asigna un factor de peso que es 
Independiente del objeto. En otras palabras si un 
impulso o fuente puntual do luz (de intensidad 
normalizada) escolocado (f.5.10) enel plano-" en 








laposición — (£:.ns) — , entonces podemos esperar 
que el plano imagen x-y responda con una función 
9%) igual an:x,;:0,,1+,Nótese que g(x,y) no es ya 


Una función puntual sino que en principio cubre todo 
el dominio de la imagen, en todo caso teniendo un 
máximo para una cierto punto (x,,;+) . Sila fuente 
Puntual deluz es ponderada por elvalor de laimagen 
deentradaer .,»,) seespera entonces quela salida 
¡comporte de la siguiente manera 








Por lo que en el límite de un número infinito de 
fuentes puntuales ¡luminandola escena (o imagen de 
entrada), se tiene 


que es la relación (5.46) del capítulo anterior, cuya 
versión digital esla ecuación (5.50) y que aquí ha sido 
obtenida — intuitivamente por consideraciones 
puramente físicas. Esto se hace aquí para 
complementar la idea que se tiene de un sistema 
linealde transformación de imágenes. Veamos ahora 
laideafísica de no-linealidad: la función de respuesta 


al impulso o función de transferencia depende de la 
¡imagen de entrada, detal manera que análogamente 
alcasolineal, supongamos que se coloca en el plano 
“1, dela imagen de entrada, en el punto (2.1) 
'unafuente puntual de luz, porlo que se espera queen 
el plano x-y de la imagen de salida se tenga una 
respuesta cuya magnitud es 


AA) 


donde (5.1) eslaimagen de entrada. Ahora bien, 
si esta respuesta es ponderada por el valor de la 
imagen de entrada en (£.1:) ., entonces puede 
esperarse que la salida, en plano imagen sea como 
sigue 





Seapreciade esto úlimo que enel planoimagen so 
tendrá toda una distribución de intensidades 
caracterizada básicamente porh, como el sistema es 
noxineal esta distribución cambiará para cada fuente 
puntual de luz, adiferencia del caso lineal.on donde la 
distribución es la misma cambiando solamente su 
valor máximo al ser ponderada por la imagen de 
entrada en cada punto de ésta. En otras palabras en 
el caso no-lineal, la función de respuesta h está ya de 
por sí modulada por la imagen de entrada. 
Imaginemos ahora que recorremos la fuente puntual 
de intensidad normalizada sobre el plano xy a 
manera de exploración pora imagen de entrada, de 
aquí que la imagen de salida sea simplemente la 
superposición de todas las distribuciones 
individuales de cada punto explorado porla fuonto de 
luz, lo que resulta entonces en la formación de la 
imagen de salida de acuerdo a la siguiente ecuación 





acon = |  10a (kan D)ALAn 





que es finalmente la ecuación (5.47) y cuya versión 
digitales muy fácil de escribi. De nuevo cuenta y de 
manera similar al caso fineal, esto complementa 
físicamente la idea que se tenía hasta ahora de 
sistemas no-Ineales. 

Enmuchos sistemas de formación de imágenes en 
dondeel espesor de la escena es de interés, setienen 
funciones de transferencia que dependen de ésta, 
siendo este el caso para radiografías y neutrografías. 
En otras situaciones experimentales puede ocurrir 
¿que haya otros objetos diferentes a la escena, con 
espesores y densidades variables interaccionando 
con ésta y que produzcan también funciones de 
transferencia no-ineales. Un ejemplo de esto último 
es la interferencia atmosférica en la captura de 
imágenes de satóite. A veces, a primera vista, un 
sistema de estos puede aparecer lineal cuando en 
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fuente 2D 










Fig. 6.15. Modelo de formación de una radiografía o neutrografía. 








realidad no lo es, establecióndose estrictamente 
¡como sistemas descritos por la ecuación de arriba; 
“aunque si bien pueden no apartarse mucho de un 
comportamiento lineal, En otras palabras, estos 
sistemas son aquellos en los cuales el fenómeno 
portador de información (ondas electromagnéticas, 
acústicas o haces corpusculares) pasa a través del 
objeto (escena) de interés, o bien cuando el 
fenómeno de transporte atraviesa un medio de 
interferencia una vez habiendo interaccionado con el 
objeto. Ejemplos de estos son los mecanismos 
formación de radiografías y la captura de imágenes 
través de la atmósfera (figs. 6.15 y8.16), Enelcaso de 
radiografías biomédicas la dependencia de la función 
de transferencia con respecto a la escena es de bajo 
orden ya que la energía delos rayos-X- empleados es 
baja; en otras palabras este sistema es noxineal pero 
sin apartarse mucho de la Iinealidad, En sistemas de 
rayos-X de ata energía como los empleados en la 
industria, donde los objetos de interés son mucho 
más densos que en biomedicina, el fenómeno de 
dispersión de Mie (Kondratyev, 1969) puede ser de 
primer orden. En el caso de imágenes obtenidas por 
satóltes orbitales, la dependencia de la función de 
transferencia con respecto a la escena puede se 
primer o segundo orden (Lira y Olwer, 1983) 
dependiendo de la longtud de onda de la radiación 
electromagnética, de la distribución de diámetros de 
las componentes atmosféricas y de las partículas en 
'semisuspensión (contaminantes) que se encuentren 
en la almósfera al momento de la captura de la 











imagen. Un posible camino de solución para 
sistemas no-lineales de segundo orden, es decir, que 
no se aparten demasiado de la linealidad, es el de 
expander en series de Taylor la función de 
transferencia alrededor an) , donde a 
es un parámetro deducible del sistema físico 
particular que se estó estudiando, entonces 








(6.175) 


Como un ejemplo de este proceso de quasi 
linealización, supóngase que la función de 
transferencia puede ser modelada como 
circularmente simátrica (gaussiana) con una varianza 
proporcional a:*(z .n jes decir, 








e pitan? 


AS A 





(6.176) 
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e imagen de salida g(x.y) 
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Fig.6.16. Modelo de formación de imagen a través de la atmóslera y por medio de un sensor satelitario. 


Mr yitam 








Desde luego, dependiendo del problema 
pañicular, otros modelos para la función de 
transferencia pueden ser escogidos. Entonces la 
primera derivada para h es 


E meso(- pre yy 


PE 





(61m 


Por lo que si se utilizan los dos primeros términos 
de la expansión de Tafcr so tencia la siguiera 
expresión aproximada para la funci 
transferencia 





e nl 
pel pr pnl, 
(6.159 

Esta última exprosión const, por tanto, de un 
término constante y uno ineal en 1 (suponiendo 


Land) + de 


NTE A 








Ea e (ng) 








podido del ración del etómero portador 
de información con el medio degradante: la oscena.o 
un medio de interferencia (figs.6-15y6.16). Amedida 
que se agregan más términos a la relación (6.177) so 
¡supone quela función de transferencia se aparte más 
de la linealidad. Al proceso arriba descrito se le 
¡conoce como quasiinealización. 
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Agentina: Noviembre 4:12, 1987. 





SEPTIMA REUNION PLENARIA SELPER, SEGUNDO SIMPO- 

(O LATINOAMERICANO DE PERCEPCIÓN REMOTA, ll Sim 
POSIO COLOMBIANO DE PERCEPCIÓN REMOTA, Coloquio 
Internacional SPOT, Primera Asunión JOBRESA (cin Board ot 
Remoto Sensing Actrles). Organizado por SELPER, nstto 
Geográfico Nacional Agustin Codazzt (IGAC) y Soctedad Colom 
lana de Fologrameta y Percepción Pamota (SCF). con. apoyo. 
¿e diversas organizaciones inamacionales (ESA. CNES, SPOT 
IMAGE, DORNIER. INPE. y otras) Bogota, Colombia Noviembre. 
1020, 1987. 








CURSO DE ENTRENAMIENTO SOBIE EL SATELITE DE MA, 
DAR EUROPEO (ERS-1), Organizado por la Agencia Espacial 
Europen (ESA) como colaboración a SELPER Baja Calomia, 
México. Noviembre 1988. 


OCTAVA REUMION PLENARIA SELPER Y TERCER SIMPOSIO 
LATINOAMERICANO DE SENSORIAMENTO REMOTO. Orgart- 
zado por SELPER yla Untwersigad Nacional Autónoma de México. 
TUNA) con el apoyo de diversas organizaciones intencionales 
fon espacial EOSAT), Acapulco, México, diciembre 5:, 1988. 


REUNION DEEXPERTOS PARA REVISION FINAL.DEL DICCIO- 
MARIO SELPER DE PENCEPCION REMOTA Patrocinada pora 


Agonci 
“gonará a SELPER esto Diccionario, producto 
“Canaaó/Poró. Pachacamac (Uma), Perú. Junio 516, 1989. 


MOVENA REUNION PLENARIA SELPER, CUARTO SIMPOSIO 
LATINOAMERICANO DE PERCEPCIÓN REMOTA, CURSO 
INTERNACIONAL SPOT Y OTRAS ACTIVIDADES ASOCIADAS. 
"Organizada por SELPER, Insttuto de Tecnología Industrial (NTD. 
Comisión Nacional de Invesigaciones Espaciales (CNI), con el 
poyo de diversas organizaciones Intemacionales. Bartoche, 
Algentina, noviembre 20:24 1989. 


OTRAS ACTIVIDADES DIVERSAS PROMOVIDAS POR LA 
SEDE SELPER/BRASIL, incluyendo co organización y copatro- 
¿nio de Cursos, Seminaros. que serán etalacas en 
las primeras publicaciones cala Rovisa SELPER en 1900. 


(S8C) ana cer natonal and interior agencies. 


INTERNATIONAL MEETING OM TRAININO NEEDS ANO USER 
ASSISTANCE Berlín, RFA, august31-soplember 6,1988, Organ- 
Led by the German Foundation lor intemational Development 
(DSE). SELPER paricpaton as inviod socaty 


"UNITED NATIONS MEETING OM SPACE SCIENCE AND TECH- 
MOLOQY IN THE EDUCATION FRAME. Co.sponsored by Mex: 
can Govemment Mxico DF. MEXICO. October 13:17, 1988, 
SELPER oia invtea. 


FIRST LATIMAMERICAN MEETINO OM SEARCH AND RESCUE! 
SATELUITES. Seo José dos Campos, SP, Braz, Octover 20-30, 
1966. Organized by INVE (Brasman Space Resoaren insfuto) and 
úcosponsored by SELPER, or the SARSAT Syntem developed by 
Canada, Francs, USA, USSA. 


“SECOND NATIONAL SEMINAR ON REMOTE SENSINA AND 
TS APPLICATIONS. Organizoo by Peruvian Soctly ol Photo- 
¡grammetry and emcto Sensing (APAFA) and co sponsored by 
SELPER Lama. Porú. May, 11-18, 1987. 


TRAININO COURSE OM EUROPEAN RADAR SATELLITE (ER: 
y. Organazoa y European Space Agency (ESA) andco-sponsored 
by SELPER Buenos Ars, Argntina. November 4-17, 1987, 








LOMBIAN SYMPOSIUM OM REMOTE SENSINO, INTERNA- 
TIONAL SPOT COLLOQUIUM, FIRST MEETINO OF JOBRESA. 
[fini Boerd ol flemote Sensing, Acimtie). Co-organizod by 
'SELPER Agustin CodazzlNaicnal Geographic insitut (GAC) and 
Colombiam Society of Photogrammetry and Remoto Sensing 
(8CP) Cosponsored by several intemalional organizations 
Bogota, Colombia. November 16-20, 1987. 


TRANNO COURSE OM EUROPEAN RADAR SATELITE (ENS: 
Y). Organized by Europesn Space Agency (ESA) lo Indorse 
SELPER actes. Daja Calfomia, Moxieo. NOvombor 1088, 


EJOMTM PLEMARY MEETINQ OF SELPER AND THIRD LATI- 
NAMERICAN SYMPOSIUM OM REMOTE SENSING, Co.organ- 
1200 by SELPER ana Nalonal Universi of Mexico (UNAM), Co- 
gonacred by several Intemaiona!organizaions. Acapulco, Mox 
Ko. Decembar5-9, 1988, 


MEETIMO OF LATINAMERICAN EXPERTISES FOR FINAL RE- 
VIEWINO OF TME REMOTE SEMSINO DICTIONARY OF 
SELPER, Sponsores by Canada Inerallonal_Davolopment 
"Agency (CIDA) ana co organized by SELPER, ONERN, CONIDA. 
Pacnacamac (Lima). Peru. June 5-16, 1989, 





ored y severa Inemational organizations. Bartocho, Argen. 
November 20:24, 1989. 


MORE IMPORTANT ACTIVITIES PROMOTED BY SELPER 
SEAT BRAZIL, including co-organization and co-sponsoring ol 
Coses, Seminars, Congress inside and cutsido Braz, ich wal 
be included in nex! SELPER Reviews (1990). 
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necesitado conocer cómo es el terreno, 
deslindes, formas, posición de edificaciones 
y accidentes nalurales, etc.; a objeto de 
planificar, ejecutar y mantener trabajos u 
obras en las distintas áreas de la vida 
cotidiana. 

La Aerofotogrametría se ha constituido en 
una herramienta eficaz, la que a travós de 
un desarrollo vertiginoso y acorde a 
nuestros tiempos, permite de manera 
fidedigna conocer la forma del relieve 
terrestre. Para esto y con la ayuda de 
modernos instrumentos. restiluidores, 
procesa las fotografías aéreas en forma 
analitica, apoyándose con la red geodésica 
nacional hasta lograr el dibujo topográfico 
del relieve a la escala de representación 
requerida, 

Sus usos, dependiendo de la escala final 
de representación, van desde 
levantamientos de predios, comunas para 
fines catastrales y extensas áreas para 
planificación y control de obras de 
ingenieria (embalses, represas, canales, 
caminos, tendido de líneas de alla tensión 
y otros). 

El Instituto Geográfico Militar posee en 
estos momentos el instrumental más 
moderno que, unido a la experiencia de 
sus profesionales, se constiluye en líder 
en este campo en Sudamérica, 


INSTITUTO GEOGRAFICO MILITAR 
Nueva Santa Isabel 1640 — Santiago-Chile 


| Desde la antigúedad, el hombre ha 


INSTITUTO GEOGRÁFICO MILITAR 
Nueva Sta, Isabel N* 1640 Fono: 6968221-8 











IMAGEN METEOSAT (36500 Km), ESA 
1492-1992; Hacia los 500 años de cooperación Europa/América 
1492-1992: Towards the 500 years of cooperation Europe/America 










IMPRESO EN LOS TALLERES DEL INSTITUTO GEOGRAFICO MILITAR DE CHILE 
Nueva Santa isabel 1640. Santiago, Chile. F: (0562) 6968221. TLX: (0562) 441677. FAX: (0562) 6088278 





